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ABSTRACT 

The RCW41 star-forming region is embedded within the Vela Molecular Ridge, hosting a massive 
stellar cluster surrounded by a conspicuous Hii region. Understanding the role of interstellar magnetic 
fields and studying the newborn stellar population is crucial to build a consistent picture of the 
physical processes acting on this kind of environment. We have carried out a detailed study of 
the interstellar polarization toward RCW41, with data from an optical and near- infrared polarimetric 
survey. Additionally, deep near-infrared images from the NTT 3.5 m telescope have been used to study 
the photometric properties of the embedded young stellar cluster, revealing several YSO's candidates. 
By using a set of pre-main sequence isochrones, a mean cluster age in the range 2.5 - 5.0 million 
years was determined, and evidence of sequential star formation were revealed. An abrupt decrease 
in R-band polarization degree is noticed toward the central ionized area, probably due to low grain 
alignment efficiency caused by the turbulent environment and/or weak intensity of magnetic fields. 
The distortion of magnetic field lines exhibit a dual behavior, with the mean orientation outside the 
area approximately following the borders of the star-forming region, and directed radially toward the 
cluster inside the ionized area, in agreement with simulations of expanding Hii regions. The spectral 
dependence of polarization allowed a meaningful determination of the total-to-selective extinction 
ratio by fittings of the Serkowski relation. Furthermore, a large rotation of polarization angle as a 
function of wavelength is detected toward several embedded stars. 

Subject headings: Galaxy: open clusters and associations: individual: RCW41 — ISM: clouds — ISM: 
magnetic fields — Stars: formation — Techniques: photometric — Techniques: 
polarimetric 
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1. INTRODUCTION 

It is well known that sites of massive star formation 
are mainly found inside heavily obscured dusty molecular 
cores, which are comm only associated to giant molecula r 



cloud complexes (e.g., Lada & Lada 2003 Lada 
Photometric studies at near-infrared (NIRj wave. 



2010[ ) 

engths 



make possible to probe deep inside where star forma- 
tion is taking place, with the very young stars frequently 
showing large infrared excesses due to the presence of 
warm dust from the associated circumstellar disks. 

On the other hand, the study of the interstellar dust 
grain's interaction with the local magnetic fields, is a 
powerful tool to understand some of the physical pro- 
cesses acting in such environments. In this context, mag- 
netic fields might play a key role in star formation. In 
fact, several physical processes, like the initial cloud col- 
lapse by the fragmentation of early stage star-forming 
cores, the channeling of interstellar material via ambipo- 
lar diffusion, the disk and bipolar outflow formation, the 
angular momentum transport, etc, possibly are strongly 
affected by the strength and configura tion of the lo 
cal magnetic field components (Mestel fc Spitzer 1956 

: de 



Nakanoil979l|Mouschovias fe Paleologou|1981[|Shu et al. 



-k Based on observations collected at the European Southern 
Observatory (La Silla, Chile), National Optical Astronomy Ob- 
servatory (CTIO, Chile) and Observatorio do Pico dos Dias, 
operated by Laboratorio Nacional do Astrofisica (LNA/MCT, 
Brazil). 
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fc Shu|[T989| IHeitsch et aL][2004| [Girart 
et al.|[2009 ). However, it is still unknown whether mag 



netic fields or interstellar turbulence provide the main 
supporting source against the cloud collapse, u ltimately 
defining th e star formation efficiency ( Padoan et al. 20041 
Crutchcr 2005) [Heiles fc Crutcher(|2O05| |McKee ^TOs 
triker,2007| ^ 

TEe formation of massive stars causes a large impact 
on the surrounding environment, ionizing the interstel- 
lar gas and providing a source of turbulence due to 
strong stellar winds and outflows from the newborn stel- 
lar population. Only recently, theoretical studies have 
attempted to describe a realistic view of the impact of 
the expansion of Hii region s on the underlying structures 
of the magn etic fle ld lines (|Kruniholz et ar]|2007 Peters 



|2011[ ). D^nfoiTunaMy 



et al.||2010[ |2011[ [Arthur eTaT 

mappings of the magnetic field structures along the ma- 
jority of Galactic Hii regions are scarcely available, al- 
though these could serve as important tests and provide 
additional constraints to such models of star formation. 

An important method to map the sky-projected lines 
of magnetic fleld using optical or NIR spectral bands, is 
from interstellar polarization by aligned dust grains, that 
serves as a well known tracer. A detailed study of linear 
polarization toward Hii regions and young stellar clusters 
can provide useful information regarding the structure of 
the Galactic magnetic fleld in these regions. 

The Vela Molecular Ridge (VMR) is an interesting 
Galactic complex where star formation is taking place 
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Fig. 1. — (a) Composite JHKs image (RGB) from the 2MASS survey showing the cluster in the direction of the IRAS 09149-4743 source, 
as well as the regions associated to our NIR photometry datasets: the dashed-line cyan rectangles indicates the regions related to Dataset 
A (J, H and Ks filters) and the solid-line cyan rectangle is related to Dataset B (H and Ks filters), (b) Images of the cluster region from 
Dataset A, at filters J {left), H {middle) and Ks {right), (c) Mosaics of the entire region of RCW41 covered by Dataset B at H {left) and 
Ks {right) filters. 
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where one can see the presence of several large interstel- 
lar features, such as the Vela supernova remnant, the 
Gum Nebula, and other structures. Beyond this "local" 
structure, at the 1 — 2 kpc range, several young star clus- 
ters are found, as well as numerous optical Hii regions 
with several sig nposts of embedded star formation such 
as H2 O masers ( Braz fc Epchteiii||1983 Zinchenko et al. 
19951 



Among these structures, a giant molecular complex ex- 
tends over 20° in the southern sky. It can be ide ntified 
by its strong CO emission, and was first studied by |Mur- 
phy & May ( |1991[ ), who divided the main structure into 
tour regions, named Aio D. Clouds towards these areas 
are located along a wide range of distances, f rom 700 to 
2000_EC pJseau et al |1992| . According to [Murphy fc" 



Maylpgi] ) regions A, C and D are located quite closer 
to us (~ Ikpc), with the cloud B being slightly more 
distant (~ 2 kpc). Recently, a s ubmillimetric survey car 



ried out by Olmi et al. (2009) detected approximately 
140 proto- and pre-steilar cores in the direction of cloud 
D. 

The IRAS 09149-4743 source is an example of young 
stellar cluster in this region. It is part of the cloud A, 
being located at an estima ted distance of 1.3 ± 0.2 kpc 
(Roman-Lopes et al. 2009|). It shows several candidate 
young stellar objects (YSOs), as well as, two massive 
stars of spectral types 09v and BOv, as revealed by previ 



ous photometric and spectroscopic infrar ed surveys ( Or- 
tiz et al. 112007 Roman-Lopes et al. 20091. Several radio 



and infrared observations are available in the direction of 
IRAS 09149-4743, tracing a number of chemical elements 



in this region, which are mainly related to the existence of 
the surrounding Hii region and the massive star forming 
area. Ortiz et al. (2007) compiled a list of all observa- 
tional data co llected from th e literature up to that date. 
Furthermore, Pirogov et al. (2007) reported the detec- 
tion of CS, N2il~'" and 1.2 mm dust continuum emissions 
toward IRAS 09149-4743. Specifically, the dust emission 
represent an almost spherical core approximately super- 
posed to the stellar cluster direction. Pirogov, ( ,2009) car- 
ried out further studies of these data through analysis of 
the cloud's radial density profile. 

Our main goal in this work is to study the stellar popu- 
lation of the associated massive star cluster, as well as, to 
investigate the structure of the interstellar magnetic field 
in the direction of the related Hii region. This was done 
by combining new deep NIR imaging with optical and 
NIR polarimetric data of the related region. The paper 
is organized as follow: Section 2 describes the NIR pho- 
tometric data and reduction techniques, as well as, the 
polarimetric observational sample. Results and analysis 
obtained from the photometric and polarimetric data are 
separately exposed respectively in Sections 3 and 4. Dis- 
cussion of these results are introduced in Section 5, and 
the final conclusions are listed in Section 6. 

2. OBSERVATIONAL DATA 

2.1. Near-Infrared Photometry 

The raw imaging data for the IRAS 09149-4743 region 
were retrieved from the ESO archive, and the observa- 
tions comprise two separate programs with the follow- 
ing identifiers: 073.D-0102 (PI Dr. Sergio Ortolani) and 
080.D-0470 (PI Dr. BenDavies). The observing missions 
were conducted during the nights of May 5th 2004 and 
February 10th 2008 using the Sofl ("Son of Isaac") NIR 
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Fig. 2. — Comparison between Sofl and 2MASS magnitudes 
as well as magnitude histograms used to derive the completeness 
limits, shown above for Dataset A (J, H and Ks filters) and below 
for Dataset B (H and Ks filters). The comparison diagrams were 
used to determine the zero point corrections for each filter, and 
the respective zero point uncertainties are indicated inside each 
diagram. The dashed lines represents the ideal correspondence 
between the two systems, and were inserted in order to clarify 
the comparison. Only isolated and good quality photometry stars 
(flag "A") from the 2MASS catalog were used. The magnitude 
histograms were used to calculate the completeness limits (cl) to 
each filter from Datasets A and B, and are indicated by the vertical 
dashed lines. 



et al. 



veyed 



camera dMoorwood et al.||1998 |), mounted on the NTT 
3.58m telescope trom La Silia/Chile. 
In Figure we use a combined 2MASS ( Skrutskie 
2006 1 image of the studied region to show the sur- 
area. The 2004 survey corresponds to the two 
small squares shown in this figure (delimited by cyan 
dashed lines); one centered on the stellar cluster and the 
other on a nearby control field. This survey will be here- 
after denominated "Dataset A", and corresponds to J, 
H and Ks bands observations using a 0.144"/pixel plate 
scale, which allows a 2.1' x 2.1' field-of-view. The obser- 
vational strategy included 24 expositions at each of the 
three filters, corresponding to 12 observations of the clus- 
ter field intercalated with 12 observations of the control 
field. Moreover, successive observations of the same field 
were made by jittering the telescope by a small displace- 
ment (w 9.5", roughly in the NE-SW direction) between 
each exposition. The FWHM values of point-like sources 



at this observational run are of ~ 1" . 

The observational data from the 2008 survey (which we 
hereafter denominate "Dataset B") was obtained using 
the H and Ks bands, with a plate scale of 0.288"/pixel, 
allowing a field of view of about 4.9' x 4.9'. Four jittering 
positions were used, moving the telescope in the N-E-S- 
W directions while positioning the cluster in each of the 
four quarters of the frame. Using this method, a larger 
field could be mapped covering an area of 7.4' x 7.4' on 
the sky, centered on the stellar cluster. In Figure [l^ the 
region associated to Dataset B corresponds to the larger 
square delimited by a solid line. Atmospheric conditions 
during this observational run allowed mean FWHM val- 
ues of « 0.7". 

Imag es from bo th datasets were treated by using 
IRAI'n [^dy|l986[ ) tasks from the NO AO package to per- 
form the following reduction steps: division by a fiat-field 
frame and illumination pattern, correction of bad pixels 
from the frame using a bad pixel mask, correction of the 
"crosstalk effect" , combination of jittered images in or- 
der to create sky fields, subtraction of the sky images, 
combination of the sky-subtracted images and alignment 
of images from different filters. The reduced images of 
Datasets A and B are respectively displayed on Figures 
[i]d and[l}:. MONTAGlQwas used to create the 7.4' x 7.4' 
field mosaic of Dataset B (Figure flfc). 

DAOFIND (from the DAOPHOT package) was used 
to locate stars peaked with 3(T above the background, 
and objects missed by this routine were further inserted 
by visual inspection. Astrometric configuration of the 
objects' world coordinate system was done by using a 
sample of 2MASS isolated objects from the same field. 
Correction of the coordinates was achieved with a rms 
uncertainty of « 0.05 arcsec in right ascencion {a) and 
declination (5), for both Datasets. 

PSF photometry was further performe d using the al - 
gorithms from the DAOPHOT package ( |Stetson|[l987[) , 
with the following photometric parameters tor each 
Dataset: PSF radius of 2.30" and fitting radius of 1.30" 
for Dataset A; PSF radius of 2.59" and fitting radius of 
0.86" for Dataset B. Several runs of psf fitting and stellar 
subtraction were performed in order to reveal and obtain 
magnitudes for the very close faint companions from the 
most crowded areas. 

Photometric calibration was achieved by comparing 
the instrumental magnitudes with a sample of isolated, 
good quality photometry stars (fiag "A") at the same 
field from the 2MASS survey. Such comparison is shown 
in Figure [2] for both datasets. These diagrams show a 
quite gooa correlation, with zero point uncertainties of 
0.02 mag in all photometric bands for both datasets. In- 
dividual magnitude uncertainties were computed using 
a quadratic sum of the instrumental errors from the psf 
photometry and the calibration uncertainty. 

Furthermore, approximate completeness limits were 
derived from the analysis of the histograms, also shown 
in Figure 2] It is defined as the point where the his- 
tograms deviate from the straight line, which repre- 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ Eletronic address: http://montage.ipac.caltech.edu/ 
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Fig. 3. — Diagrams related to Datasets A (a) and B (b), used to verify a possible color dependence associated to the photometric 
calibration of the SofI magnitudes based on the 2MASS magnitudes. The horizontal dashed lines were inserted for clarity. The dispersion 
shows that no presence of color dependences are evident within our instrumental uncertainties. 



sent a linear fit to the logarithmic of the number of 
objects per magnitude bin. The values obtained for 
Datasets A and B are (J,H,Ks)Hmzt = (18.6,17.4,16.8) 
and (R,Ks)iimit = (18.4,17.8), respectively. Consider- 
ing the above photometric limits of both datasets, an 
improvement of 1.6, 2.1 and 2.8 magnitudes have been 
obtained respectively in J, H and Ks bands, as compared 
to the earlier survey by Ortiz et al. (2007), allowing us 
to probe deeper into the lower-mass stellar population. 

In order to check the existence of color dependence 
terms to be applied to the SOFI photometry, we con- 
structed diagrams of magnitude differences (between the 
SofI and 2MASS systems) as a function of the (J-H) and 
(H-Ks) colors, as can be seen on Figures [3^ (Dataset A) 
andlsb (Dataset B, only the H-Ks color). As can be no- 
ticed7 there is no need to apply color correction terms to 
the derived SOFI photometry. 

The photometric results are shown in the Appendix, 
respectively corresponding to 530 and 2608 stars in 
Datasets A (Table [l]) and B (Table |2|. The full versions 
of the tables are available only in tne electronic version 
of this paper. 

2.2. Optical/ Near- Infrared Polarimetry 

The optical and NIR polarimetric data used in this 
work have been separately obtained with different tele- 
scopes and instrumental configurations. The optical data 
were collected using the 0.9 m telescope of the CTIO 
(Cerro Tololo, Chile, operated by the National Opti- 
cal Astronomy Observatory - NOAO) observatory, dur- 
ing February and April/2010; and March/2011. The 
NIR data were obtained using the 1.6 m telescope of the 
OPD (Pico do Dias, Brazil, operated by Laboratorio Na- 
tional de Astrofi'sica - LNA/MCT) observatory, during 



April/2011. 

The polarimetric modules used in both telescopes are 
similar and consist basically of a rotatable achromatic 
half-wave retarder followed by a calcite Savart plate and 
a filter w heel (for a complete d escription of this instru- 
ment, see Magalhaes et al.|1996 |. The half- wave retarder 
can be rotated in steps of 22.5, and one polarization 
modulation cycle is covered for every 90° rotation of this 
wave-plate. This arrangement provides two images of 
each object on the detector, which correspond to the 
perpendicular polarizations beams fo (ordinary) and /e 
(extraordinary). By rotating the half- wave plate by 45° 
yields in a rotation of the polarization direction of 90°. 
Thus, at the detector area where fo was first detected, 
now /e is imaged and vice versa. Combining all four 
intensities reduces flat-field irregularities. In addition, 
the simultaneous imaging of the two beams allows to ob- 
serve under non-photometric conditions and, at the same 
time, virtually suppress the sky polarization component. 
Besides, if extended interstellar emission from the field 
contributes with some polarization level, this component 
is also automatically eliminated. Finally, each set of a 
polarimetric observation was performed using blocks of 
8 steps of the half-wave retarder, which corresponds to a 
180° rotation of the plate. 

Reductions were performed in the standard manner us- 
ing IRAF's routines, for both the optical and NIR data. 
Point-like sources were selected from the images using 
DAOFIND, for stars peaked 5cr above the local back- 
ground level, with the saturated objects removed from 
the sample. The total counts for each source were com- 
puted applying aperture photometry using the PHOT 
routine for several ring sizes around each star. Objects 
affected by cosmic rays, bad pixels or superposition of 
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Fig. 4. — The image above is a RGB combination of the R (DSS), 
J and Ks (2MASS) bands, showing both the Ha extended emission 
related to the Hll Region (blue) and the position of the embedded 
stellar cluster. The areas mapped by the polarization survey are 
indicated by the boxes related to different spectral bands. The R 
optical band (yellow solid line) has been used to cover a large area 
(28' X 28') encompassing both the cluster and the Hll region, while 
V and I optical bands (green dashed box) cover a 8' X 8' area more 
focused on the cluster. The NIR observations with the H band (red 
dashed-dotted box) cover a 3.5' X 7.0' area towards the cluster and 
extended along the south. 

light beams from close companions were not used in our 
analysis (classified as bad polarization values). 

From the difference in the measured flux for each beam, 
we derived the polarimetric parameters using a set of 
IRAF tasks specific ally designed for this purpose (PC- 
CDPACK package, [Pereyra| 2000| . This set includes a 
FORTRAN routine that reads the data files and cal- 
culates the normalized linear polarization from a least- 
square fit solution, which yields the degree of linear po- 
larization (P), the polarization position angle (9, mea- 
sured from north to east), and the Stokes parameters Q 
and U, as well as the theoretical (i.e., the photon noise) 
and measured errors. The latter are obtained from the 
residuals of the observations at each wave plate position 
angle (ipi) with respect to the expected cos4'0i curve. In 
further analysis of the polarimetric data, we adopted the 
greater value between both estimated errors. 

Finally, in order to determine the reference direction 
of the polarizer, and to check for any possible intrinsic 
instrumental polarization, a set of polarimetric standard 
stars were observed at each night. The obtained polariza- 
tions degree for the observed unpolarized standard stars 
proved that the instrumental polarization are negligible 
for both telescopes and instrumentation. 

In the case of the optical survey, V, R and I Johnson- 
Cousins' filters (covering different areas over the studied 
region) were used. In Figure El it is shown a false-color 
image measuring about 30' x 30' that corresponds to a 
combination of the R (DSS=blue), J(2MASS=green) and 
Ks (2MASS=red) images of the RCW41 region. The ex- 
tended emission seen in blue is probably mostly due to 



the Ha emission of the associated Hii region; the clus- 
ter objects are mainly seen in red and near the centre. 
The polarimetric R-band survey (yellow box - 28' x 28') 
covers almost the whole area, while the V- and I-bands 
covers about 8' x 8' of the region centered on the cluster 
(green dashed line). Each field was observed using expo- 
sures of 300 sec per half- wave plate position, and in order 
to cover the entire area of the R-band survey, a 3 x 3 
frames mosaic mapping was conducted at the CTIO's 
0.9 m telescope. The polarized standard stars for the op- 
tical survey were selected from the catal ogs compiled by 
Tapia (1988) and |Turnshek et"ar (19901, corresponding 
toTElD 2US^3, HD 111579, HD 1265STand HD 110984. 

The NIR observations gathered at the 1.6 m telescope 
of the OPD observatory were performed using the H- 
band of the NIR camera CamlV, which has a 4' x 4' field- 
of-view when mounted at this telescope. At each position 
of the half-wave plate, sixty 10 sec images were obtained 
(in order to avoid the detector's non-linearity limit for 
the brighter sources), jittering the pointing by 12 arcsec 
between each image in a cross-pattern, summing up to a 
600 sec of total exposure time. Two fields were mapped 
in this manner, one centered on the cluster and a second 
one shifted towards de south, therefore covering a 3.5' x 
7.0' area as shown in Figure [4j by the red dashed-dotted 
box. The NIR polarimetric standards Eli as 14, Elias 25 



(p Oph ) and HD164740 were selec ted from Whittet et al 
|l992| , [Wiiking et al.| (|l980l), and [Wilkin g et al.| ( |1982| ).' 

I'he complete polarimetric survey is presented in the 
Appendix (TablelSl) with only a portion of the table given 
here just to demonstrate its format and content. The 
complete Table is given only in the electronic version of 
the article. 

3. RESULTS AND ANALYSIS FROM THE DEEP NIR 
PHOTOMETRY 

3.1. Color-Color Diagrams 

Valuable information regarding the nature of the clus- 
ter's stellar population may be inferred from (J — H) x 
{H — Ks) diagrams, which are shown in Figure p] They 
were constructed from Dataset A, for objects in the clus- 
ter and control areas detected in all three bands. There 
are represented the unreddened main sequence (MS), 
along with the giant/supergiant branch locus (shown at 
the left bottom b y the red solid line), as obtained from 
Koornne"e^ fl983) and corrected to the 2MASS p hotomet- 
ric system using the transformation proposed by|Carpen-| 



ter 



( 2001 ). Also, the locu s of the Classical T Tauri (CTT) 



stars (|Meyer et al."|T 997[) is indicated in the cluster's dia- 
Dv the solid da 



gram by the solid dark green line, wit h the standard red- 
dening law (Rieke & Lebofsky 1985[ ) being represented 
(in both cluster and control diagrams) by the red dashed 
lines. Different colors and symbols are used to represent 
objects from different parts of the diagram. 

By comparing the distribution of sources in both dia- 
grams, we can notice that despite the majority of control 
field stars show extinction levels similar to the cluster 
field stars, the objects from the cluster region present a 
much larger number of highly reddened objects and stars 
with color excess (respectively red asterisks and orange 
triangles). Furthermore, by analyzing the distribution of 
the majority of points from the cluster's diagram (mainly 
the gray open circles), we note that they are approxi- 
mately distributed along a strip (indicated by the solid 
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Fig. 5. — Color-color diagrams for the cluster (a) and control 
(6) areas from Dataset A. Stars marked with Objl-6 label s are ob - 
jects previou sly studied through spectroscopic techniques (Roman- 
|Lopes et al.|[2009| ). The red dashed lines represent the recldenihg 
band trom tne standard reddening law, and the dark green lines 
denote the locus of CTT Stars. Different symbols and colors desig- 
nate different sections from the diagram (a): the bulk of cluster's 
stars is represented by gray open circles, with 1.1 < J — H < 2.2, 
within the reddening band; also within the reddening band are 
stars presenting higher (red asterisks, J — H > 2.2) and lower 
(pink diamonds, J — H < 1.1) extinction values; purple plus signs 
and orange triangles respectively designate objects towards the left 
and right (therefore with color excess) of the reddening band. A 
small dispersion of H — Ks = 0.1 beyond the limits of the redden- 
ing band is allowed, to account for photometric uncertainties. All 
objects from the control area (6) are marked with green crosses. 



blue line polygon) which is roughly parallel to the CTT 
locus (the dark green line). On the other hand, stars 
from the control field are more "vertically" distributed 
in the diagram, resembling the reddened Main Sequence 
locus. This may indicate that most of the stars from the 
cluster region are probably Pre-Main Sequence (PMS) 
stars. 

3.2. Mean Visual Extinction 

It is possible to take advantage of the fact that, the 
bulk of stellar objects from the cluster field lay in a dis- 
tribution which is roughly parallel to the CTT line locus 
(blue polygon in Figure pa) , to estimate the mean visual 
extinction in the direction of the cluster. We initially 
will assume that the objects located inside the polygon 
area are CTT stars. Therefo re, assuming the standar d 
interstellar reddening law of Rieke & Lebofsky (1985), 
we have applied the de-reddemng vector computing the 
visual extinction for each star inside the polygon. The 
average of these values yields a mean visual extinction of 
Av = 7.6 ± 2.0. We represent in Figure [5^ the reddened 
CTT locus {Av = 7.6), as well as the uncertainty from 
this computation, i.e., the Av = 9.6/5.6 CTT loci (dark 
green dot-dashed lines). Although a large scattering of 
the CTT band is seen, we estimate that about 30% of 
all stars from the cluster field are located within these 
limits. Besides, the interstellar extinction is probably 
spatially non-uniform along the cluster's region and con- 
tamination from foreground objects is inevitable, which 
contributes to increase the observed scattering. 

An independent computation of Ay may be inferred 
from spectrosco pically confirmed /B stars. According 
to the survey by 
as Objla, Obj2 



Roman-Lopes et al. (2009), stars labeled 
and Ubj4 in i' lgure jS^, probably belong 
to the group of cluster's most massive objects (respec- 
tively BOv, 09v and B7v). Assuming that they have 
already reached the MS, we can obtain another value 
for the cluster's mean visual extinction by de-reddening 
their observed (J-H ) and (H-K s ) colo rs, using the intrin- 
sic colors given by Koornneef ( 1983) and the interstel - 
lar extinction law taken from 'Rieke & Lebofsky (1985). 
Adopting this method, we have found the following Ay 
values, respectively for Objla, Obj2, and Obj4: 7.8±1.4, 
7.0 ± 1.8, and 6.4 ± 2.3. All values are in the range 6- 
8 magnitudes, which is consistent with the computation 
obtained from the CTT candidates method. 

3.3. Spatial Distribution of the Cluster's Stellar 
Population 

Considering that the photometric properties of the 
cluster's area indicate that it is composed by stellar pop- 
ulations with different characteristics (embedded PMS 
stars, foreground and background objects, stars with 
high extinction, etc), it is important to study their spa- 
tial distribution, relative to the cluster's location. Figure 
[6]shows a false color image of the region obtained by com- 
bining the J/H/Ks images of Dataset A. Each detected 
star is marked with colors and symbols according to its 
location on the color-color diagram, as defined in Figure 
[5^, with the stars designed as Objl-6 indicated by the 
green labels. Some very faint stars from the field were 
not marked because they were not detected in all three 
bands, and therefore were not assigned a location in the 
color-color diagram. 



Optical/NIR Photo-polarimetric study of RCW41 



7 



RCW41 Cluster Region - Dataset A: 
J/H/Ks RGB Combination (SofI Images) 




36' 3D' 




44' 42' 
Ri^hl Ascensiim (2000) 

Fig. 6. — Combined RGB image of the cluster region using J 
(blue), H (green) and Ks (red) images from Sofl's Dataset A. Con- 
trast has been enhanced in order to reveal the bright extended emis- 
sion surrounding the embedded stellar cluster. Stars are marked 
with colors and symbols related to their location on the color-color 
diagram (see Figure[5^), and objects indicated as O bi 1-6 were pre- 
viously studied through spectroscopic techniques l |Roman-Lopes| 
[et a l. 2009). The white-colored contours are from the MINUU survey 
by l^inchenko et al.| | |2dd0| , representing the 10o,io — 9o,9 transition 
witn tne intensity levels starting from 15% of the peak intensity 
(1.5 K- km s~^), raising in steps of 7.5%. 

The cluster's region is composed by two distinct struc- 
tures: the main portion, which is associated to an ex- 
tended emission visible mainly in the Ks band, is formed 
by a number of objects concentrated around the Objla/b 
stars (these two central objects present an angular sepa- 
ration of only 1.3 arcsec and are displayed in a close-up 
view at the top of Figure |6]); the small sub-cluster region 
is located 1.1 arcmin towards the SE of the main cluster's 



center, and hosts t he Obj2, which is a 09v star (Roman 
Lopes et al. 2009), presu mably the most ma ssive object 



m the held. According to Ortiz et al. (2007), Objla and 



Obj2 are the best candidates to ionization sources of the 
associated Hii region. 

The contours shown in the map (colored in white) are 
fro m the HNCO su r vey o f massive Galactic dense cores 
by Zinchenko et al. (2000), which covers mainly the cen- 
tral parts ot the cluster. The HNCO molecule is most 
easily excited radiatively (rather then collisionally) inside 
the molecular clouds' densest regions {n > lO^cm"'^), 
being particularly sensitive to the far-infrared radiation 
field from warm dust, and therefore is a suitable tracer of 
regions at the vicinity of massive star- formi ng sites. As 
previously pointed out by Ortiz et al. (2007), the HNCO 
survey revealed an area supposedly swept out by stellar 
radiation around Objla/b, as well as the presence of a 
high-density region 30-50 arcsec NW of the central stars 
(Objla/b). 

We may correlate this information with the location of 
stars with high extinction (marked with red asterisks). 



These highly obscured objects present Ay values in the 
range 15 - 35, and are most commonly found toward the 
North of the field shown in Figure [61 From a total of 
19 high extinction objects, 16 (i.e., 84% of them) are lo- 
cated above the cyan-colored dashed line which divides 
Figure [6] in Northern and Southern halves, suggesting 
an extinction rising gradient in the South-North direc- 
tion. Furthermore, 9 of these stars (representing 47% 
of them) are located within the HNCO contours, indi- 
cating the expected correlation between the high den- 
sity medium and the presence of highly obscured stars. 
Therefore, although background stars may be present, 
some of these are most probably highly embedded young 
cluster's sources, and deserve special attention. 

There is no clear trending for objects marked with gray 
open circles or orange triangles, being scattered along the 
entire area. However, for pink diamonds (i.e., lower ex- 
tinction stars or also massive members of the cluster), 
there is a weak evidence that an anti-correlation exists 
relative to the positions of high extinction stars (red as- 
terisks): 63% of the pink diamonds (17 among 27 stars) 
arc located within the Southern area of the image. Al- 
though it is not possible to assert this unequivocally, this 
fact corroborates the idea that a rising extinction gradi- 
ent occurs in the South-North direction. This evidence 
may support the idea that interstellar material surround- 
ing Objla is only beginning to be swept out (according to 
the structure of HNCO contours), while the area around 
the subcluster (which hosts the 09v star) have already 
been cleaned due to the action of its strong winds. 

3.4. Color-Magnitude Diagrams 

In Figure [T] we show three color-magnitude diagrams 
for both the cluster and control areas. While the con- 
trol field is probably composed by foreground and back- 
ground main sequence and giant stars, the cluster's field 
also contains members of the cluster itself. In order to 
locate the overlap region, stars from the control area 
(denoted by green crosses) were plotted in all diagrams. 
This is useful to identify stars from the cluster area pre- 
senting photometric characteristics of contaminating ob- 
jects. 

Figures ^{Hx [H-Ks] diagram) a.nd^ {Ksx[J~Ks] 
diagram) were both created with colors and magnitudes 
from Dataset A, u sing the cluster and control areas as 
did in Section 3.1 (the inset in Figure [7^ shows the clus- 
ter's and control areas overploted in a small 2MASS Ks 
band image). Points from these two diagrams are plot- 
ted using the same color/symbol scheme as in Figure 
[5] therefore allowing us to identify the positions of the 
different groups previously defined in the color-color dia- 
gram. The zero age main-sequence (for solar metall icity) 
is represented by a 2MASS JHKs Padova isochrone ( Bon- 
atto et al.|2004 1 , which is denoted by the black solid line 



surrounded by a gray band that accounts for the ±2.0 
mag uncertainty in Ay- As a complement, we also plot- 
ted three pre-main-sequence (PMS) isochrones {t = 1, 
2.5, and 5 Myr, for solar metal licity) taken f rom young 



stellar evolution models by ^SiessetaLl ( 2000 1 . We used 
the cluster's distance {d = 1.3 ± 0.2 kpc) and the in- 



terstellar extinction as in Section |3.2[ Furthermore, an 
arrow corresponding to an amount of 5 mag in visual 
extinction, was also included in each diagram. 
Figure [fj; is a iJ x (iJ — Ks) diagram built with the H 
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a) H X (R-Ks) Diagram (Dataset A) - Cluster field stars: O >K O A + ; Control field stars: x 




O * 

'^^^ .^1 



b> Ks X (J-Ks) Diagram (Dataset A) - Cluster field stars:0 J« O A + ; Control field stars: x 



a 14- 




FlG. 7. — Color-Magnitude diagrams using Dataset A (a: H X [H — Ks]; and b: Ks X [J — Ks]) and Dataset B (c: H X [H — Ks]). The 
small 2MASS images inside diagrams a and c show the cluster and control areas used for each Dataset. Stars from the control areas are 
marked with green crosses, while cluster's stars in diagrams a and b are marked with colors and symbols according to its location on the 
color-color diagram (Figure [51. All stars from the cluster area in diagram c are denoted by gray open circles. The MS and several PMS 
isochrones are indicated, ana the yellow line provides a statistical separation between the contaminating field stars and objects from the 
cluster, as shown in Figure [S] In dia gram b, the inset image shows the location of the cluster's objects from the turn-on area (blue arrows), 
i.e., located within the blue dashed rectangle. Obj2 was saturated in Dataset B, and therefore is not shown in diagram c. 
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and Ks data from Dataset B, which provides a larger area 
around the cluster (see Figure [T]). As shown by the inset, 
we have kept the same clusters area used in Dataset A 
(Figures [7^ and b) and chosen a different control region, 
encompassing a wide area at the external borders of the 
region covered by Dataset B (defined by the green cross- 
hatched). This should allow a more complete probe of 
the field stellar population. Since J band data are not 
available from Dataset B, a position in the color-color 
diagram was not assigned, and therefore, we have simply 
used gray open circles to represent cluster's stars and 
green crosses for objects from the control area. 

By analyzing the distribution of the cluster's and con- 
trol's stars in all three diagrams, we notice some striking 
common features. Initially, the majority of the point 
sources from the cluster's region appear shifted from the 
MS locus, being displaced to its right with higher color 
values, i.e, over the region where the PMS isochrones are 
present. This indicates that the majority of the cluster's 
sources are probably very young stars that still have not 
reached the hydrogen burning phase. However, as would 
be expected, some objects from the cluster area are lo- 
cated toward the left part of the diagram, occupying the 
same area as the control field stars, and therefore are 
probably members of the Galactic field stellar popula- 
tion. Also, from figures Wk and b, we can see that ob- 
jects positioned toward the left side of the ZAMS are 
mainly those represented by pink diamonds and purple 
plus signs, i.e., lower extinction stars and objects with 
anomalous colors, therefore indicating that such objects 
are likely foreground stars. On the other hand, it is ev- 
ident that some of the pink diamond brighter stars (for 
example. Objects la, 2 and 4), are in fact cluster's mem- 
bers that have already reached the MS phase. 

Although the bulk of the cluster's stars (represented 
by the circles), are located amongst the PMS isochrones, 
a certain number of objects appear distributed to the 
right in the diagrams, therefore presenting large J — Ks 
and H — Ks values. This is probably a combined ef- 
fect of intrinsic infrared excess from PMS stars, non- 
uniform interstellar extinction along the clusters' area, 
and photometric uncertainties. We specifically notice 
those stars represented by the orange triangles, which 
are well displaced toward the right showing high color 
indexes, which is an indication of emission from circun- 
stellar disks. Furthermore, these objects are relatively 
faint with H > 16, therefore probably corresponding to 
the very low-mass component of the YSO candidate sam- 
ple. 

Finally, the objects represented by red asterisks (iden- 
tified as objects within the reddening band, but show- 
ing high extinction levels) appear presenting some of the 
highest J — Ks and H — Ks values. Although some of 
these stars may be background sources, considering the 
non-uniform and clumped nature of the interstellar envi- 
ronment around the cluster (see, for example, the HNCO 
emis sion contours in Figure m\ and the discussion in Sec- 



3.5. The Cluster's Evolutionary Status 



tion |3.3 ), it is possible that some of them are actually 
highly embedded cluster's PMS stars. For example, it is 
highly probable that Objlb (the close companion to the 
central star Objla, and that is marked with a red aster- 
isk) is actually a massive YSO. More on this subject will 
be further discussed in Section l3^ 



The PMS isochrones from Siess et al. ( 2000 ) shown in 
Figure [7] may be correlated with the stellar distribution 
in the color-magnitude diagrams, to provide a mean es- 
timate of the cluster's age. However, contaminating field 
stars from the cluster's area should be previously identi- 
fied, at least in a statistical manner. Such identification 
is based on the overlap areas on the color-magnitude di- 
agrams between the cluster and control populations. 

In order to provide a more quantitative recognition of 
the contaminating objects, the following method have 
been adopted: the three color-magnitude diagrams from 
Figure [7| were divided in bands of 1 magnitude, and 
to every strip, a histogram was constructed, using the 
color index of each respective diagram. The results are 
shown in Figure |8] where histograms from parts a, b, 
and c are respectively related to the color-magnitude di- 
agrams from Figures [7^, ^)p, and [Tb. Each histogram is 
divided between the cluster (red solid lines) and control 
(blue dashed lines) populations, and bins of 0.15 mag 
were used. The histograms were normalized by the total 
number of stars used in the color-magnitude diagrams 
of each area, namely: A^tot (Dataset A, cluster) — 302, 
A'^tot (Dataset A, control) = 202, A^tot (Dataset B, cluster) 
= 384 and Atot (Dataset B, control) — 396. Therefore, 
these diagrams represent the percentage fraction of stars 
within each magnitude strip and color bin. 

A rough separation between the intrinsic cluster's stars 
(mainly displaced toward the right of the histograms) 
and the field stellar population may be defined by the 
color index value where the cluster's fractional number 
of stars becomes predominant relative to the control. To 
every histogram presenting a sufficient number of points, 
this value have been marked with a yellow vertical line. 
Therefore, it represents a separation criterion between 
the field stars (objects from the left of the yellow line) and 
the intrinsic cluster's population (stars from the right of 
the yellow line). These values were also indicated in Fig- 
ure [7] by the yellow connected dots and lines. This pro- 
vides an appropriate scheme to roughly de-contaminate 
the color magnitude diagrams, therefore allowing a more 
clear analysis of the isochrones' positions relative to the 
distribution of possible intrinsic points from the cluster. 

Specially when analysing Figures |7^ and[7j: (the H x 
[H — Ks] diagrams), we notice that tne yellow separation 
lines passes exactly between the locus related to t = 2.5 
Myr and t — 5.0 Myr. Moreover, in these two diagrams, 
the bulk of the cluster's points (gray open circles), which 
are toward the right of the yellow line, seems to roughly 
define a band around the t — 2.5 Myr isochrone. 

On the other hand, the turn-on point is better defined 
in the Ks x {J — Ks) diagram (Figure [7]d) , as indicated 
by the blue dashed rectangle, which encompassed both 
objects that have already reached the MS, and those that 
are currently leaving the PMS stage. Notice that over the 
MS region, below the turn-on point for the t = 5.0 Myr 
isochrone, there is a small gap for cluster's sources with 
magnitudes in the range from Ks — 13 up to Ks ~ 14 
(from whereon contaminating objects seem to dominate), 
suggesting that there are no stars from the cluster with 
t > 5.0 Myr. Furthermore, the Obj4 star, whose dis- 
tance is consistent with be ing an intrinsic cluster's source 
( Roman-Lopes et al.|2009 ), is located over the MS, at the 
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a) Dataset A: H strips, H-Ks histograms 
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b) Dataset A: Ks strips, J-Ks histograms 
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c) Dataset B: H strips, H-Ks histograms 
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Fig, 8. — Histograms of color indexes for cluster (red solid lines) and control areas (blue dashed lines), separated by successive strips of 1 
magnitude wide. Histograms from parts a, b and c are directly related to the distribution of the stellar population in the color-magnitude 
diagrams respectively from Figures [7^ '[Zt'' ^iicl[7|;. Therefore, parts a and b are both from Dataset A, respectively representing histograms 
oi H — Ks and J — Ks in magnitude divisions or H and Ks. Histograms from part c is from Dataset B, where the H — Ks color is used 
in bands of H magnitudes. Each value in the histograms were normalized by the total number of stars used in each field. Vertical yellow 
lines are used in several histograms to help identifying contaminating stars from the cluster's region, by separating the overlapping areas 
between control and cluster populations (at the left of the yellow line) from the population which is mainly composed by cluster's stars (at 
the right). 
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midway between the t = 2.5 and t — 5.0 Myr isochrones. 

These evidence suggest that the cluster's mean age is 
between 2.5 and 5.0 Myr. Moreover, these analysis re- 
veal that a sequential star formation scenario may have 
occurred within this region. The spatial location of ob- 
jects from the turn-on area in Figure [7j3 (within the blue 
dashed rectangle) is indicated by blue arrows at the in- 
set image in the same diagram. It is evident that most 
of these objects (10 among a total of 13) are positioned 
closer to (or within) the subcluster, instead on the central 
main cluster area (as would be expected if mass segre- 
gation processes had already taken place). Therefore, a 
plausible explanation to this evidence may be that the 
most massive objects (located at the subcluster), were 
formed earlier, subsequently inducing star formation at 
the main cluster area. 

This evolutionary view is supported by the fact that 
the interstellar material near the subcluster seems to 
have already been swept out, while at the clus ter area 
denser clumps are still present (see Section 3.3). If the 
multi-epoch formation hypothesis is in fact true, the large 
spread of cluster's points on the color-magnitude dia- 
grams towards higher color values could also be explained 
by an intrinsic age spread of the studied stellar popula- 
tion. 

3.6. Analysis of Individual Objects with known 
Spectroscopic Properties 

The spectroscopic f eatures of several spe cific objects 
previously studied by Roman-Lopes et al. ( [2009| ), may 
now be correlated with the objects' position on the color- 
color and color-magnitude diagrams. The same stars 
from that survey were renamed as Objl-6 and are in- 
dicated in Figures [5} |6j and [7) 

Three of the most massive objects from the cluster - 
Objla, 2, and 4 - were respectively classified as BOv, 
09v, and B7v-B8v stars. Their location on the color- 
magnitude diagrams (Figure [t]) are consistent with stars 
that have already reached the MS. 

Objlb sho ws spectral fines ch aracteristic of YSO's, and 
according to Ortiz et al. ( |2007 ), the set Objla-|-b also in- 
dicates features ot war ni dust, such as infrared emission 
beyond 5^m (probably coming from Objlb). In fact, this 
object is located within the high-extinction region of the 
color-color diagram, an indication of its embedded na- 
ture. Furthermore, we notice that its K -band spectrum 
(see Figure 6 from Roman-Lopes et al. ( 2009j )) presents 
features which are very similar to the spectrum of T Tau 
Sb, the thir d member from the T Tau system (Duchene 



et al. 2002): both sources show a large Br7 emission 
(which IS a strong accretion signature), as well as several 
absorption CO overtones. This evidence suggests that 
Objlb, the close companion to the massive BOv central 
object (Objla), is a very young T Tauri star. The pres- 
ence of this young object near the center of the main 
cluster region provides additional support to the idea 
that star formation within this area have been triggered 
by the earlier star for matio n activity at the subcluster, 
as discussed in Section 1331 

Although the sources labeled as Obj5 and Obj6 were 
previously classified as late-type field stars, theirposi- 
tion on the color- magnitude diagrams (Figures wk, h, 
and c) are consistent with a classification as T Tauri 
stars. Even though their spectroscopic features may be 



interpreted as from typical low mass late-type stars, ob- 
servations have been reported of Class I stars with few 
spectroscopic signs of youth. Two examples are the IRAS 
03220-1-3035(8) and IRAS F0325 8+3105 sources, from 
the N IR spectroscopic survey by |Connelley fc Greene] 
(2010), which were classified as young YSO's, despite 



presenting K-band spectra only with some CO overtones 
in absorption and almost no Br7 emission, hence being 
very similar to Obj5 and Obj6. Therefore, if we also take 
into account their position at the very crowded area near 
the cluster's center (see Figure [6]), these sources are most 
likely T Tauri stars, members of the cluster. 

Other interesting properties are related to Obj3, that 
is a highly reddened source, and has NIR colors corre- 
sponding to a visual extinction of Ay > 20, which is con- 
sistent with its spatial location among the HNCO con- 
tours (Figure [6]). Although the previous analysis of its 
K-band spectrum showed that it could be a typical late- 
type background star, its position in the color-magnitude 
diagram suggests that it may also be interpreted as a 
highly reddened medium mass YSO. A careful analysis 
of its spectrum show a weak set of CO absorption lines, 
which is usually a feature of low-mass YSO's. However, 
its high NIR luminosity rather suggests a more massive 
nature. Therefore, a valid explanation for the weakness 
of these lines could be due to veiling of the photosphere 
and CO b ands by the circunstella r dust emission, as sug- 
gested by [Casali fc Eiroa| (jl995|). In this case, an anti- 
correlation of these lines' intensity with the NIR color 
excess is expected. Although more observations are nec- 
essary to fully characterize this source, these evidence 
indicate that Obj3 is likely a medium mass YSO from 
the cluster. 



Bin— 5 arcsec 




-47.88 
-47.90 
-47.92 
-47.94 
-47.96 
-47.98 
-48.00 




9.285 9.282 9.279 9.276 9.272 
Right Ascension (hours) 



Bin= 15 arcsec 



9.285 9.282 9.279 9.276 9.272 
Right Ascension (hours) 



Bin= 20 arcsec 




-47.88 [lRA..Dcc.xq27 3 33h,47 94018 dcfi) 
'ivi 

-47.90 \ 



-48.00 




I 



9285 9.282 9.279 9.276 9 272 
Right Ascension (hoLirs) 



9.285 9.282 9.279 9.275 9272 
Right Ascension (hours) 



Fig. 9. — Stellar density distributions to different bin sizes, cre- 
ated using the photometric positions of Dataset B. The gray-scale is 
related to the stellar density, where darker tones represent a higher 
density, as can be inferred by the scale bars. Bi-dimensional Gaus- 
sian fits applied to each distribution are represented the red contour 
lines, which are related to the following levels; 0.013, 0.015, and 
0.030 stars/arcsec^. The center positions computed from the Gaus- 
sian fits are specified at the top of each diagram, and are indicated 
as a red bullet. 
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Fig. 10. — (a) King profile fitting using concentric rings around 
the cluster's center position. The values obtained for cro, Rcore, 
<^ fields a-nd Riim are expressed at the top of the diagram, and the 
spatial size of Rcore and Rum are indicated as dashed red circles 
in the Ks image from Dataset B (inset). The error bars represent 
the Poisson uncertainties. The afield parameter, which represents 
the uniform stellar density from the field, is separately computed, 
by taking the mean value (and the standard deviation) of points 
with r > Riim- (b) The same King profile fitting, ignoring the blue 
cross-hatched area, which contaminates the RDP due to a stellar 
density scattering toward the South-East. Notice the sharper for- 
mat of the curve fitting in this case. In both diagrams (a) and 
(b), the thick yellow line represent the radial density profile of the 
1.2 mm dust continuum emission (250 GHz) originated from the 
cluster's dire ction, defining an almost spherical pro jec ted core, as 
described by |Pirogov et al.| (20071 and |Pirogov] i |2009| l . The scale 
at the right represent tne 25U Litiz intensity, and tne peak of this 
radial density profile is adjusted to match the peak of the stellar 
density from the King function, for comparison purposes. 

3.7. Cluster's Center Determination and Radial 
Density Profile 

In order to perform a rigorous analysis of the cluster's 
Radial Density Profile (RDP), it is important to deter- 
mine its mean center. Therefore, we have used the Ks 
image from Dataset B, divided the region into bins of 
specific widths (5, 10, 15, and 20 arcsec), and computed 
the stellar density distribution. The result is shown in 



Figure [9j where the gray-scale represents the number of 
stars per bin area. 

The cluster's center position for each bin width was 
determined through bi-diniensional Gaussian fits to the 
stellar density distributions, represented by the red 
contour maps superposed to each diagram in Figure 

fA mean value has been computed, resulting in 
^^,4)=(9hl6'°43.5%-47°56'24.2"), with an uncertainty 
of (0.2'', 0.9"). Although the cluster's stars appear 
slightly scattered toward the South and South-East (near 
the subcluster), a highly peaked concentration appears 
around UcS^, the cluster's center. 

The morphology of the RCW41 cluster's RDP may 
be s tudied by fitting s of a two-parameter King func- 
tion ( King|1962 1966), which may be approximately de- 
scribed in terms of the central density {(Jq) and the core 
radius (Rcore)- This was achieved by positioning several 
circumcentric rings around the computed center, with 
widths of 10 arcsec each, and computing the stellar den- 
sity as a function of the distance from the center {a{r)). 
The obtained density profile was fitted by the following 
equation: 



cr(r) = a field + 



1 + {r/Rcoref 



(1) 



The results are expressed in Figure [Tok, where the 
points indicate the stellar densities at each ring, and the 
King function fitting is denoted by the solid line, together 
with the uncertainty curves (dashed lines) obtained from 
the uncertainties of the fitted parameters. Starting from 
the center, the stellar density decreases towards larger 
radius, and the point where the cluster's density com- 
pletely merges with the field density is visually identified 
as Riim, i.e., an approximate value for the cluster's lim- 
iting radius. 

Note that a bump extending from r w 65 to 110" ap- 
pears in this radial profile, which is due to the stellar 
scattering towards the south-east area (which includes 
the sub-cluster - small blue box in Figure flOk). In ¥ig- 
ure |10| 3, we have removed this area from the RDP anal- 
ysis (blue cross-hatched region), in order to provide a 
cleaner fitting only of the cluster's main concentration, 
i.e., the stellar densities were computed from increasing 
half-rings at the Northwest area of the field. The new 
King function is much sharply concentrated toward the 
center in this case, showing a large decrease of 44% in 



and 38% in i?, 



lini • 



The estimated value of R, 



lini 



Re 

indicates that the control areas used in the photometric 
analysis from the previous Sections, are sufficiently dis- 
tant from the cluster's main structure, therefore avoiding 
si gnificant contaminat ion. 

|Lada fc Lada (2003) pointed out that embedded stellar 
clusters may be morphologically classified as hierarchical- 
type (with multiple peaks at the density profile, fre- 
quently found scattered over a large spatial scale) or cen- 
trally condensed- type (highly concentrated distributions, 
with relatively smooth RDPs). Although the density dis- 
tributions show that the cluster do present some isolated 
stellar concentrations, the relatively smooth fit of the 
King function indicates that it is probably of the cen- 
trally condensed-type. This structural type is a common 
feature of very young stellar clusters, since most of the 
stars have not yet dispersed throughout the field, and are 
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still heavily concentrated around the center. 

Such morphology may be correlated with 250GHz ob- 
servations of the 1.2mm dust continu um emission to- 
wards this region (Pirogov et al. 2007), that defines an 
almost spherical core m which the cluster is probably im- 
mersed in. Based on these observations, Pirogov| ( 2009 1 
computed a power law radial density prohie ot the dust 
emission, which is shown in Figures 10 1 and \10}p by 
the thick yellow line (the peak emission was adjusted to 
match the maximum value of the King function). One 
may note that both distributions (from the dust emission 
and from the stellar density) are similar in size, suggest- 
ing that the stellar system is still embedded within part 
of the original cloud in which star formation was initially 
ignited. Therefore, the presence of such dense interstel- 
lar material increases the gravitational potential and so 
the physical bounding of the system, aiding to keep the 
stars spatially concentrated. The centrally condensed- 
type stellar cluster occurs because of the interstellar gas 
and dust from the parental molecular cloud, which are 
still present in the cluster's region, and probably begin- 
ning to be swept away due to the action of the new-born 
massive stars. 

It is important to point out that in the above men- 
tioned comparison, the centers of both the dust emission 
radial profile and the King function, are not exactly the 
same, being separate d by about ~ 11 arcsec. Further- 
more 
this o 



Pirogov ( 2009 ) adopted a much larger distance to 
3ject (2.6 kpc), as compared to the value of 1.3 
kpc sp ectroscopically determined by |Roman-Lopes et al.| 
(2009). However, the difference is irrelevant in this case, 



since we have converted the spatial dimensions to angular 
sizes in Figure [TO] 

4. RESULTS AND ANALYSIS FROM THE 
OPTICAL/NEAR-INFRARED POLARIMETRY 

4.1. The Large-scale Distribution of R-band 
Polarization Vectors 

The R-band polarization survey covers a large area 
around the RCW41 region, encompassing both the stel- 
lar cluster and the ionized area, as previously shown in 
Figure |4j A polarization degree and angle were assigned 
to every point-like source detected in the field, and by 
this manner we have constructed a map of the Hii region 
showing the distribution of polarizatio n v ectors super- 
posed to each star, as exposed in Figure [TT] The vectors' 
orientations indicate the predominant direction of the 
electric field oscillations related to the stellar light beam, 
and assuming the standard grain alignment mechanisms, 
the vectors trace the sky-proj ected component of the in - 
terstellar magnetic field lines ( Davis &: Greenstein|195i ). 

In this image, the size of each vector is proportional 
to the polarization degree (a 10% sized vector is indi- 
cated at the top-right of the graph as reference), and the 
vector's colors are related to the quality of the polar- 
ization measurement: good quality measurements, with 
P/ap > 3, are shown in yellow, while lower quality data, 
with 3 > P/ap > 1, are shown in red. Green circles 
denote objects with P/ap < 1 or stars with bad polar- 
ization values (i.e., affected by cosmic rays, bad pixels of 
the detector or light beam superposition with close com- 
panions). Therefore, all optically detected objects from 
this survey are marked, either with a vector or a green 
circle, according to the above mentioned criteria. 



Since the observations were done by integrating the 
light of all stars from the field with the same expo- 
sure time, the classification of an individual measure- 
ment among one of these groups depends on the stellar 
brightness (which increases the individual signal-to-noise 
ratio) and also on the intrinsic polarization degree from 
the source (therefore, unpolarized sources naturally have 
larger probability to fall on the P/ap < 1 classification). 

Considering that many stars are too obscured to be de- 
tected by the optical survey, only the brightest sources 
were visible. Even though, these observations are still 
useful to study the large-scale properties of the magnetic 
field structure around the Hii region, as will be shown 
in the following sections. However, it is important to 
point out that, considering the 1.3 kpc distance to the 
cluster, a fraction of these stars are foreground objects, 
and therefore their observed polarization levels are possi- 
bly mainly caused by interstellar components not related 
to the RCW41 region. These properties will be further 



analyzed in Section 4.2 



The background image used to produce Figure 11 
shows interstellar structures related to the star-forming 
region extending into a large area around the embed- 
ded stellar cluster. The image is a RGB combination of: 
the DSS R-band image (blue), showing optically visible 
stars and the ionized region (the Ha emission is within 
the spectral region defined by the R band); the 2MASS 
Ks image (green), indicating the location of several ob- 
scured stars visible only in the near-infrared, including 
the associated stellar cluster, nearly at the image's cen- 
ter; the MSX 8.2^m mid-infrared emission (red), showing 
mainly the radiation from the warm dust and polycyclic 
aromatic hydrocarbons (PAHs). 

A striking feature seen in this image is that the mid- 
infrared emission (red) defines a great ring or bubble of 
warm dust and PAHs that encompasses the ionized area 
denoted by the extended emission in blue. This structure 
suggests that UV radiation from massive stars and high 
energy photo-ionized electrons from inside the Hii region 
heats the dust grains and excites the PAHs from the sur- 
rounding environment, generating the "envelope" type 
of emission. It is not uncommon to find mid-infrared 
emission of PAHs at the outer edges of Hii regions, as 
revealed for exampl e by the Spitzer GLIMPSE su rvey of 



the Galactic Plane ( Churchwell et al. 



2006 2007) 



Some mixture between the R-band and mid-infrared 
emission is seen in the direction of the cluster, where 
several interstellar filaments extend almost radially from 
it. Such morphological structure will be further analyzed 
in Section |431 

Another interesting information is provided by the 
IRAS 100/im emission which is attributed to the colder 
dust component and is represented in Figure [TT] by the 
gray contours. Note that the 100/im emission shows 
a peak almost centered on the cluster's area and that 
the outer contours (especially the 170 MJy/sr isophote) 
roughly follow the morphology of the hot dust ring from 
the mid-infrared radiation, suggesting the existence of a 
cold dust "cocoon" around the entire area. 

Comparing the global distribution of polarization vec- 
tors with the direction of the Galactic Plane (the b = 0.8° 
line), we note that both directions are not correlated. 
All-sky optical polarization surveys have shown that, 
particularly along the Galactic Plane, polarization vec- 
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RCW41: R band vectors, R(SDSSyKs(2MASS)/8.2 ^{MSX)RGB image 

P/{iP>3 ■ 3>P/{TP>1 C P/ of < 1 & bad P values 

IRAS 100 contours, levels: 128, 170, 240, 5(K) and l(X)OMJ>/sr 
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Right Ascension (2O00) 



Fig. 11. — Distribution of R-band polarization vectors at the RCW41 region. The image is a 30 X 30 arcmin RGB combination of the 
DSS R-band (blue), the 2MASS Ks band (green) and MSX 8.2/im emission (red). Notice the large mid-IR ring or bubble (red) surrounding 
the blue-colored ionized area. The embedded cluster is located nearly at the center of the region. Vectors' sizes are proportional do 
the polarization degrees and colors are related the data quality, as specified above the diagram. Green circles denote objects with poor 
{P/cTp < 1) or bad polarization values. Gray lines are the IRAS 100/.tm contours, showing the 128, 170, 240, 500, and 1000 MJy/sr emission 
levels. 



tors tend to present an overall horizontal orientation, i.e 



fc Ell is 
noted 



By 



1976 



Heiles 

latter 



parallel to the plane (^Mathewson fc Ford 1970 Axon 



Crutchcr 
authors, 



2UU5p. However, as 
)oTari 



lose this tendency on the magnetic "poles", located at 
I « (80°, 260°). The VMR is located along one of such 
"pole" (the RCW41 cluster is at Z = 270.3°), and there- 
fore the correlation between polarization vectors and the 
Galactic Plane direction is indeed not expected. 

The distribution of polarization vectors presents some 
characteristics that can be related with the features of 



the Hii emission region. Note that toward the ionized 
area, and specially near the cluster, the vectors' sizes are 
small when compared to the outer vectors. This evidence 
indicates a decrease of polarization degree values toward 
the center of the region. Furthermore, at the right side 
of the field, a global bending of the vectors' directions is 
seen both at the area's bottom-right and top-right por- 
tions, resembling to roughly follow the curvature of the 
IRAS lOO/im emission contours in these regions. Such 
properties will be further analyzed in Section [4. 3[ 
It is important do point out that, since a fraction of the 
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Fig. 12. — Histograms related to the R-band polarimetric dataset, 
representing the polarization angle (top) and polarization degree 
(bottom) distributions. Yellow and red colors respectively denote 
the poor and good-quality data, while the black histograms con- 
sider the entire sample presenting P/ap > 1. 

sample is composed by YSO's with infrared excess, some 
of these sources may contribute with a contaminating 
polarized emission from the circunstellar disk. However, 
the large-scale correlated pattern of polarization angles 
seen in the map may not be explained by intrinsic polar- 
ization. The overall characteristics are best described by 
interstellar polarization due to dust dichroic absorption, 
with some possible scatter in polarization angles due to 
individual disk emission from some sources. 

Figure [12] shows histograms of polarization angle and 
polarization degree divided by different P/dp intervals. 
The Gaussian fits to the angle histograms show that 
the dispersion (indicated by the FWHM lines) of the 
3 > F/cTp > 1 data (red) is larger than the P/ap > 3 
case (yellow), an expected trend given the lower data 
quality. However, both distributions are quite well corre- 
lated, showing polarization angles predominantly in the 
160° — 10° range. Moreover, by comparing the orienta- 
tions of higher and lower quality vectors in Figure |11[ 
it may be noted that these are quite well correlated lo^ 
cally, indicating that the lower quality data may also be 
included in the analysis. Polarization degrees varies be- 
tween and ~ 16%, with a large number of stars peaking 
at « 2 - 6%. 

4.2. Correlation with 2MASS photometry 

Since the individual stellar distances for the polarimet- 
ric survey are not available, we shall take advantage of 
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Fig. 13. — (a) (J — H) x {H — Ks) color-color diagram us- 
ing 2MASS data for the R-band polarimetric sample; (b) Pr{%) 
x(J — H) diagram constructed from a correlation between the R- 
band polarization data and the 2MASS magnitudes; (c) Spatial 
distribution of probable foreground (green crosses) and background 
or embedded stars (orange vectors, represented as asterisks in dia- 
grams (a) and (b)). 



the 2MASS data to approximately infer the foreground, 
background and intrinsic cluster's populations. Given 
that the Hii region is located at 1.3 kpc from the sun, it 
is important to evaluate if the interstellar polarization of 
objects from such distance range are significantly affected 
by the nearby structures along the same line-of-sight. 
Therefore, the goal is to identify, at least statistically. 
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foreground objects and analyze their polarization levels 
in compar ison with the typical values from the field. 

Figures 13 a, and b show respectively the (J — H) x 
{H — Ks) color-color and the Pr x (J — H) diagrams us- 
ing the R-band polarimetric sample and the associated 
2MASS data. Figure [13^ shows that the bulk of polari- 
metric points present \ j—H) values in the range 0.5—1.1 
mag, in contrast with the previous photometric analysis 
(Figurejsj), where it has been shown that stars both from 
the cluster and control regions are somewhat more red- 
dened. This is probably due to a combination of two 
factors: (1) the earlier analysis was carried out toward 
a much smaller area, concentrated on the cluster, where 
the Ay values are statistically higher; and (2) the polari- 
metric survey is optically limited, and therefore only the 
brightest sou rces were detected. 

Figure [13) 3 shows that, for those objects having J — 
H > 0.9, there is almost no star presenting Pji < 1%. 
This effect is emphasized by the vertical {J — H — 0.9) 
and horizontal (Pr = 1%) green dashed lines in this dia- 
gram. This is a suggestion that objects obscured by the 
molecular cloud (and therefore located behind or within 
it) show a minimum polarization level equal to 1%. As a 
consequence, stars presenting polarization values below 
this limit are probably foreground objects. Such points 
are marked with green crosses, and define a thick band 
around the unreddened MS locus in the color-color di- 



agram (Figure 13 1) , a further evidence that they are 
mainly foreground, not absorbed stars. Moreo ver, their 
spatial location on the field is uniform (Figure 13:), as 
expected for a randomly distributed group of foreground 
objects along the line-of-sight. 

We conclude that although a small polarization level 
due to foreground structures is probably included in the 
stellar light for most of the objects in the observed field, 
this level is small when compared to the typical polariza- 
tion levels found for the region (which is mainly between 
2 and 6%, as suggested by th e polarization distrib ution 
given in Figure 12 1. In fact, Santos et al. (2011) have 
shown that in theairections near / ~ 260" over the galac- 
tic plane, polarization degrees in the V band are smaller 
than 1% at least up to a distance of 400 pc from the sun. 
Therefore, it is expected that a major fraction of the po- 
larization vectors for the observed field indeed map the 
magnetic field lines within the Hii region. 

By analyzing Figure[l3^, we note that a group of points 
(marked with orange asterisks) are located above a small 
gap (at J — H K, 1.1) within the reddening band of the 
color-color distribution. The Ay = 7.6 MS reddened 



locus (as calculated for the cluster. Section 3.2 1 is super 



posed to these points, showing that these stars are red- 
dened by the molecular cloud, and therefore are probably 
either background or embedded objects. Their spatial 
distribution is also shown in Figure |13[ :. Open circles, 
which comprises the main group of stars from the field, 
have typical polarization degrees above the 1% level for 
foreground objects and typical J — H values below the 
1.1 estimated limit for background and embedded stars. 
Therefore, this group is probably composed by a mix- 
ture of foreground, background, and mainly embedded 
objects in the Hii region. 



4.3. Large-scale Mean Polarization Degree and Angle 

Maps 

In order to analyze further the evidence that the degree 
of polarization is lower inside the Hii region, we have cho- 
sen an approximate center of the ionized area (based on 
the R-band extended emission, represented by the blue 
color in the image shown in Figure fTTl), and studied the 
behavior of the mean polarization as a lunction of this ra- 
dial distance from the center. This analysis is displayed 
in Figure |14^ : beginning from the selected center posi- 
tion ((a, dj^ooo = (9 16 42.9, -47 59 48)), several circum- 
centric rings of increasing radius were built, and inside 
each ring, the Q and U Stokes parameters from individ- 
ual objects (with P/ap > 1) were averaged, therefore 
allowing to compute the mean polarization value associ- 
ated to that ring and its associated standard deviation. 
Three different ring widths were used in this analysis 
(0.75, 1.00 and 1.25 arcmin), and the results are identi- 
fied with distinct symbols in Figure [14^. Therefore, each 
mean polarization value was plotted as a function of the 
angular distance from the chosen center. 

A sharp increase in mean polarization degree (from 
« 2.8% to sa 3.9%) was found at a distance of ~ 7.8 
arcmin from the center of the Hli region (correspond- 
ing to ^ 2.9pc, assuming that the cluster is located at 
1.3kpc from the sun). Such distance is denot ed by an 
orange circle in the RCW41 image of Figure [T4^ . By 
comparing the position where the polarization degree in- 
creases (with respect to values inside the Hii region) with 
the MSX 8.2/im mid-infrared emission (red), we find a 
striking evidence that such position roughly corresponds 
to the warm dust and PAHs "bubble" that encloses the 
ionized gas. 

Given these evidences, we have attempted to build a 
"mean polarization image" of the RCW41 region, based 
on the R-band survey. To achieve this result, we have 
divided the spatial region where the polarimetric data is 
distributed into 50 x 50 intervals of equal width. Consid- 
ering each observation inside a specific division, a mean 
polarization value was computed (only objects having 
P/(jp > 1 were used). In order to account for cases 
where no stars were available inside a particular division, 
a smoothing routine was applied, returning the same im- 
age, but averaged locally by using a small box that runs 
through the entire arrajPl The result is shown in Figure 



[T4p , where the polarization image is superposed on the 
actual picture of the RCW41 Hii region. Redder colors 
correspond to lower polarization levels while bluer colors 
indicate higher polarization levels, as specified by the col- 
orbar given at the top of the diagram. The location of 
the cluster (and the associated sub-cluster), as well as 
of the ring where mean polarization is found to sharply 
increase, are also indicated. 

Note that nearby the cluster (and slightly toward the 
northeast of it) , we find an area where polarization mean 
values are lower, in comparison with typical values for the 
field, presenting P < 2% (marked in red). Furthermore, 
toward the left and bottom side of the area (southeast), 
we find a large region where the polarization degree is 
higher {P ^ 6%, marked in dark blue). This region 
marks the transition border between the inner and outer 
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Fig. 14. — (a) Analysis of the radial dependence of mean polarization degree inside the RCW41 ionized area, where the center has been 
chosen based on the R-band extended emission (as shown by the blue color at the image toward the right). An abrupt rise in polarization 
degree is found roughly coinciding with the position of the warm dust/PAHs ring seen at the MSX 8.2/im mid-infrared emission (red), (b) 
Mean polarization map, built by averaging the polarization values of point-like sources into an image of 50 X 50 equal divisions. Notice 
that mean polarization nearby the cluster are smaller than the typical values from the field, and a large area at the bottom left following 
the ridge of the Hll region shows an increase of polarization values, (c) Mean polarization angle image, following the same procedure used 
in part (b). A global bending of the sky-projected magnetic field lines is found at the bottom right of the image, also following the rim of 
the ionized area. 



sides of the southeast part of the Hii region, and is prob- 
ably responsible for the sharp increase in polarization 
detected in Figure 14 1. 

Two models may be proposed to explain the decrease 
of polarization toward the ionized region: (1) given that 
temperature values of the local interstellar medium and 
the radiation field from hot stars are probably higher in- 
side the Hii region, it is expected that mechanisms of 
grain alignment are disturbed by the turbulent medium, 
hindering the alignment efficiency and therefore provid- 
ing lower P values. Polarization efficiency may also be 



affected if the intensity of the magnetic field inside the 
photo-ionized area is intrinsically weaker, which may be 
true if field lines have been pushed aside by the expansion 
of the Hii region; (2) dust grains have been swept away 
from inside the Hii region, due to the action of intense 
stellar winds from the newborn stars, therefore resulting 
in a lower dust column density and consequently lower 
P levels. Obviously, if the second model is true, a simi- 
lar decrease of interstellar extinction as a function of the 
distance from the center would be expected. Further dis- 
cussion on the distinction between both models will be 
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carried out in Section O 

We have also constructed a "mean polarization angle 
image" , using the same procedures described above. The 
result is shown in Figure |14| :, where yellow and red col- 
ors represent vectors slanted to the right and blue colors 
are regions with mean polarization angles slanted to the 
left (with respect to the vertical orientation). By this 
diagram, we are able to highlight global changes in po- 
larization angle with respect to the typical values found 
in the field. In fact, we note that toward the bottom 
right of the image (southwest) , a large area marked with 
green/yellow/red colors indicate a global slanting of the 
polarization vectors toward the right. Furthermore, such 
regions seems to roughly follow the outer con tour of the 
ionized area, indicated by the circle of Figure 14 1. Such 
inclinatio n o f the polarization vectors may al so b e seen 
in Figure [TTj as pointed out earlier in Section |4?T] 

This structure suggests that the physical processes that 
generated the Hll region, followed by the expansion of 
the ionized area, caused a global deflection of the orig- 
inal magnetic field lines in the field. Such deflection of 
the magnetic field lines around the expanding ionization 
front is most prominent toward the southwest area sur- 
rounding the Hii region, as exposed in Figure [14} :. As 



a) 



an example, Matthews et al. (2002) describes another 



region were a similar distortion ot magnetic field lines 
was observed toward NGC 2024, with polarization vec- 
tors being s wept due to the ex pansion of the associated 



Hii region. Tang et al 



effect toward G5.89, a 



(2009) also described a similar 
U.Ol pc sized ultra-compact Hii 
region at a much earlier evolutionary stage (^ 600 years) , 
already displaying significant distortion of the magnetic 
field lines as a consequence of the expansion of a shell-like 
structure. 

4.4. Wavelength Dependence of Polarization and 
Determination of the Total-to-Selective Extinction 
Ratio 

It is possible to take advantage of the polarimetric sur- 
vey in multiple spectral bands to analyze specific proper- 
ties of the dust particles along this region. For instance, 
the grain size distribution causes a direct influence on the 
underlying interstellar law, which is generally assumed to 
follow a standard behavior. Therefore such analysis will 
enable us to directly test this hypothesis, at least in the 
direction of the cluster and nearby surroundings, were V, 
R, I and H polarization data overlap (see Figure [4]). 

It is well known that an empirical law, known as the 
Serkowski relation, may be used to describe the spectral 
dependence of the interstellar polarization ( j Serkowski 
etaL][T975l ): 



Pa = P. 



, cxp 



-if In 



An 



A 



(2) 



where the typical value of the parameter K is 1.15, and 
Pmax and \max dcuotc respectively the maximum polar- 
ization value and the wavelength where such maximum 
point is reached. Furthermore, it is also known that when 
stellar light trespasses different dust layers with different 
grain alignment conditions, it is possible that the result- 
ing beam will present a spectral dependence of the polar- 
ization angle, us uallv represented as a rota t ion of with 
respect to A"^ ([Gehrels fc Silvester||l965[ |Coyne||T974 
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Fig. 15. — (a): Example of the adjusted Serkowski relation 
(P X A^^) and the wavelength dependence of polarization angle 
id X A~^) for a particular object presenting V, R, I and H data. 
The values of Pmax, ^max, Rv ) S'lid for this specific case are 
indicated at the top of the diagrams, (b): [P/ Pmax) x (Xmax/^) 
diagram, along with the Serkowski curve using points of all avail- 
able bands for each of the 34 fitted objects. Stars with VRIH or 
VRI data are distinguished respectively by blue or red open circles. 



Messinger et al.||1997 ). 

By analyzing the overlap region between the polari- 
metric surveys in the V, R, I and H bands, we have 
selected stars with VRIH or at least VRI observations, 
where the requirement P/crp > 3 was fulfilled in all ob- 
served bands. Our polarimetric sample contains 34 ob- 
jects meeting these conditions. Values of Pmax, ^max, 
and a measurement of the rotation of the polarization 
angle were obtained for each of these stars. As an exam- 
ple. Figure [15^ illustrates the results obtained for one of 
them, i.e., tEe P x A^^ graph (along with the adjusted 
Serkowski curve) and the 9 x A~^ diagram, also showing 
the fitted values for this particular case. 

Figure[l5]3 shows a diagram of (P/Pmax) x (A.max/A), 
where the mdividual Pmax and Xmax fitted values were 
used together with the polarimetric data for the available 
bands. Red circles indicate polarimetric data to which 
the fitting was performed only with the VRI bands, while 
blue circles include also the H band. Note that, compar- 
ing the points with the position of the Serkowski curve, 
there is a larger dispersion toward lower A values (i.e., 
toward the V band) . This is an expected effect since the 
signal-to-noise ratio of the polarimetric measurement for 
a particular embedded object is supposed to be smaller 
toward bluer colors. However, the overall fitted points 
forms a band around the Serkowski curve, therefore re- 
vealing a quite good adjustment. 

We estimated the ratio of the total-to-selective extinc- 
tion, Ry, by applying th e well known relation between 
this parameter and Xmax ( Serkowski et al.|1975 Whittet 
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Fig. 16. — (a): Distribution of the fitted values of the total-to- 
selective extinction ratio (R), along with a Gaussian fit indicating 
a mean R value of ij = 2.96 it 0.42. (b): Spatial distribution of the 
stars studied v^ith the multi-band polarization data (represented by 
the colored symbols), along with the R-band polarization vectors, 
both overlaid in a R/J/Ks combined image of the RCW41 region. 
Symbols denote different angle rotation (A9) intervals, and differ- 
ent colors represent different strips of the R value, as specified by 
the caption above the image, (c): Diagram of the obtained R val- 
ues as a function of the NIR colors, showing no particular trend or 
concentration toward lower or higher values. 



fc van Breda||1978|[Whittet|[2003| , 



Rv = (5.6±0.3)A,; 



(3) 



Figure [TBp, shows a histogram of the obtained Ry val- 
ues, indicating that most of them are concentrated be- 
tween ^ 2.5 and ^ 3.5. By adjusting a Gaussian fit to 
this distribution, the peak value suggests a mean ratio 
of the total-to-selective extinction of i? = 2.96 ± 0.42, 



which agrees with the standard value of 3.09 for the typ- 
ical interstellar grains. This finding gives support to the 
distance of 1.3 ± 0.2 kpc for the embedded ste llar clus- 
ter, obtained th rough spectroscopic techniques (Roman- 
Lopes et al.|2009). Furthermore, it also corroborates the 



position ot the PMS isochrones and main sequence locus 
in the color-mag nitud e di agra ms from the photometric 
analysis (Sections 3.4 and 3.5), since the standard inter- 
stellar la w w as assumed in those cases. 

Figure [16]3 shows a R(DSS)/J/Ks(2MASS) composed 
image of the RCW41 region, with the local R-band po- 
larization vectors distribution, and colored symbols indi- 
cating the objects from the multi-band analysis. Dif- 
ferent colors for the symbols are related to different 
strips of Rv values (divided in intervals of 0.5, as spec- 
ified above the image), and the symbols indicate dif- 
ferent intervals of the polarization angle rotation (A0, 
where □ represent a small rotation level - between —5 
and 5°//im^^ - while x and o denote high rotation 
levels respectively of -40 < A6l(°/^m~^) < -20 and 
20 < A6l(%m-i) < 40). 

Note that there is no spatial concentration of specific 
Ry or A6 values over the observed area. Although most 
of the A6 values show low rotation levels, a great A6 
dispersion is found, with extreme high-rotation values of 
—31 and 37°/^m~^. Along these directions, stellar light 
has probably passed through more than one interstellar 
layer with different grain alignment conditions. The A9 
dispersion is probably a consequence of a turbulent and 
fragmented medium associated to the Hii region, with 
superposing filaments depending on the chosen line-of- 
sight. 

Figure [TBb shows that there are no particular trend of 
the fittealtv values to any NIR color index. If stars 
with color excess due to circunstellar disks also show an 
intrinsic polarization level, this should interfere with the 
fits of the Serkowski relation, and therefore with the ob- 
tained Rv values. However, these diagrams show no bias 
of the Ry values toward higher color values. 

4.5. The Small-scale Distribution of H-band 
Polarization Vectors 

The H-band polarimetric survey is focused on the 
RCW41 cluster (Figure |4]), and since most of the objects 
in this area are highly obscured in the visual wavelengths, 
the NIR survey provides a much larger number of polar- 
ization vectors at least in a region of a few parsecs around 
the star-forming locus. 

Figure [17] shows the distribution of the H band po- 
larization vectors surrounding the stellar cluster, super- 
posed on a RGB combination of the H and Ks images 
(respectively blue and green, from Sofl's Dataset B), and 
the 5.8 /im image (red) collected using the Infrared Ar- 
ray Camera (IRAC)|^ onboard Spitzer. The morphology 
of the mid-infrared extended emission is dominated by 
several conspicuous filaments emanating almost radially 
from the cluster. These structures resemble the pattern 
previously noted in the MSX 8.2 /im emission (Figure 
11 ), but obviously with a much higher resolution, allow- 
mg the distinction of individual filaments near the clus- 
ter. 

* The data has been retrieved from the NASA/IPAC Infrared 
Science Archive (IRSA). 
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Fig. 17. — H band polarization vectors map of an area of approximately 3.5' X 7.0' encompassing the RCW41 stellar cluster. The image is 
a RGB combination of the H and Ks images from the Sofl's Dataset B NIR survey (respectively blue and green), and the 5.8^m image from 
Spitzer's IRAC (red), showing several filaments bended toward the cluster. Blue vectors, yellow vectors and green open circles respectively 
represent data with P/ap > 3 (good quality), 3 > P/crp > 1 (lower quality) and P/ap < 1 (also indicating bad polarization data). Stars 
marked with open purple squares are possible contaminating field objects, presenting the following photometric features: H — Ks < 0.5 
and H > 13.5. Histograms denote the distribution of polarization angles for objects within each of the white dotted areas labeled as 1-6 
(blue and yellow lines from the histograms also denote respectively good and lower quality data). The peak of each Gaussian fit is defined 
as the mean polarization orientation from that specific field. This mean angle is used to draw a grid inside each selected area, therefore 
denoting the mean local polarization orientation along that line-of-sight. 



Polarization vectors are plotted with different col- 
ors according to the data quality (blue vectors have 
P/(7p > 3 while yellow vectors represent data with 
3 > P/ap > 1), and green circles denote lower qual- 
ity polarization measurements. We have also adopted 
in this case a criteria to roughly separate contaminating 
field stars from cluster's members, based on the fact that 
the polarimetric survey in the H band reaches an approx- 
imate photometric limit aX H ^ 16. We have found this 
information by correlating the polarimetric data with the 
photometric Dataset B and analyzing the H magnitudes 
distribution for the polarimetric sample. Therefore, by 
studying the color-magnitude diagram of Figure [7|:, we 
notice that up to iJ < 16, the contaminating stars are 
dominant \i H — Ks < 0.5 and H > 13.5. Stars with 
such photometric characteristics are marked with a pur- 
ple square in Figure |17[ indicating that these are pos- 
sibly field stars. A consistent trend in this diagram is 
that the number of field stars within the cluster is small, 



compared to the outer regions of the image, where this 
number increases. 

The first striking evidence found in the distribution of 
polarimetric vectors is that the polarization angle dis- 
persion is much larger than in the R band large-scale 
survey. This is probably an indication of the existence of 
a largely turbulent medium near the cluster. In fact, if 
we analyse the histogram of polarization angles for the 
entire H band sample (Figure 



18 



top panel), we can 
see that the entire range of possible angles is covered. 
However, a peak is reached at similar angles as in the 
large-scale distribution, i.e., at 6* ~ 160° suggesting that 
even in smaller scales, there is an underlying trend to 
follow the global scale pattern. Furthermore, no ma- 
jor differences are noted between the distributions of the 
good quality {P/ap > 3, blue) and the low quality data 
(3 > P/ap > 1, yellow). The polarization degree distri- 



bution (Figure 18 bottom panel) shows that the major- 
ity of stellar objects present Ph values within 1 — 4%. 
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Fig. 18. — Distributions of polarization angle (top) and polar- 
ization degree (bottom) for the NIR polarimetric survey in the H 
band. Blue lines are related to the data with P/crp > 3 and yellow 
lines represent data with 3 > P/crp > 1. 

Even when considering the large angle dispersion, we 
find th at along even smaller regions within the area of 
Figure 17 some local trends may be unveiled, specially 
toward the interstellar filaments seen in 5.8/xm extended 
emission. In order to reveal these underlying features, 
we have selected six sub- areas, mainly encompassing the 
most prominent filaments, as indicated by the white dot- 
ted quadrangles marked with the 1-6 labels. Considering 
the polarization vectors within each of these areas, the 
associated histograms of polarization angles were built, 
being that objects marked as possible contaminating field 
stars were removed from the analysis. The distributions 
are also shown in Figure |17[ linked to each respective 
area, and the division between good quality {P/up > 3) 
and low quality data (3 > P/ap > 1) are indicated re- 
spectively by the blue and yellow histograms. Gaussian 
fits were applied to each distribution, and the peak of 
the adjusted function corresponds to the main trend of 
the polarization angles within that area. By using this 
main value, a white dotted grid was drawn within each 
area, therefore corresponding to the local predominant 
polarimetric orientation. 

We note that along the majority of the chosen areas, 
the main polarization direction seems to be consistent 
with the direction of the associated filament. Such trend 
is specially evident toward areas 2, 3 and 4. Specifically 
toward area 2, although only eight vectors contribute to 



the histogram, all of them are oriented parallel to the 
filament, presenting polarization angles around « 110°. 
Such orientation is almost perpendicular to the main 
trend for the large-scale pattern, indicating a clear corre- 
lation with the local filament structure. Toward areas 1 
and 6, a small difference seems to exist between the main 
polarization angle and the filament direction, although 
the orientation pointing in the direction of the cluster is 
evident. The existence of substructure or superposed fil- 
aments with different orientations along the line-of-sight 
may contribute to this displacement. Over area 5, which 
encompasses both the subcluster and a part of the main 
cluster region, an orientation pointing toward the cluster 
is also evident. 

We conclude that despite the large angle dispersion, 
magnetic field lines on a small scale seem to be defiected 
inwards, pointing toward the cluster and probably con- 
nected to the direction of the large 5.8/im filaments from 
the Spitzer survey. The polarization angle dispersion is 
indeed expected to exist in such a dynamical star forming 
region, where strong stellar winds from massive objects 
and possibly outflows from newborn lower mass stars, 
should play an important role as sources of turbulence to 
the surrounding interstellar medium. As a consequence, 
grain alignment efficiency and magnetic field lines are 
supposed to be affected, therefore characterizing a higher 
dispersion of the polarization orientations. 

5. DISCUSSION 

5.1. Comparison of the magnetic field structure with 
theoretical studies on Hll regions' expansion 



Simulations conducted by Peters et al. (2010) describes 
the collapse of a massive molecular cloud followed by the 
formation of a small stellar cluster, which ionizes the ac- 
cretion flow and causes the Hii region to expand. |Spitzer| 
( 1978 1 first described the evolution of a photo-ionized re- 



gion around an ionization source in a uniform medium, 
and provided an analytical solution to the expansion of 
the Hii region. However, several physical mechanisms, 
such as heating and cooling, turbulence, the effects of 
ionizing and non-ionizing ultraviolet radiation. X-rays, 
and magnetic fields, must be considered in order to pro- 
vi de a realistic view of the expansion process. 

Arthur et al.l (2011) have carried out simulations of 
this type considering the effects of a previously magne- 
tized medium. Comparison of these simulations on non- 
magnetized and magnetized environments provide simi- 
lar results in terms of the general morphology of the Hii 
regions, therefore suggesting that the presence of mag- 
netic field lines does not cause a significant impact in 
providing resistance to the expansion of the ionized area, 
although this influence is most evident when considering 
the morphology of small-scale features, such as globules 
and interstellar fllaments. 

On the other hand, such expansion is expected to cause 
a large impact on the overall shape of the magnetic field 
lines. The main result from these simulations have shown 
that after the expansion of the Hii region, a dual behavior 
of the field's morphology may be noticed when compar- 
ing the inside and outside regions of the ionized area: 
outside the Hii region (i.e., in the location of the neutral 
gas, following its borders), magnetic field lines tend to 
lie along the ionization front, forming a ring around the 
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ionized area; inside the Hii region (i.e., within the ionized 
area), magnetic field fines are expected to fie perpendic- 
ular to the ionization front, therefore pointing towards 
the ionization source. 

These qualitative predictions are consistent with the 
results of the large and small-scale polarimetri c su r veys 
of t he R CW41 region introduced in Sections |4.1[ |4.3[ 
and |4.5| The R-band survey have revealed that mag- 



netic field lines are bended along the borders of the Hii 
region toward some areas (Figures [Tl| and 14), which is 
consistent with the interpretation oimagnetic field lines 
being pushed aside due to the expansion of the region, 
therefore lying along the ionization front on the neutral 
areas. 

Moreover, our H-band polarimetric survey showed 
that, although a large angle dispersion is evident over 
the entire area, the mean polarization direction resem- 
bles a radially oriented field toward the cluster, roughly 
following the direction of several filaments from the mid- 
infrared Spitzer image. This is consistent with the theo- 
retical prediction that inside the ionized area, field lines 
are expected to lie perpendicularly to the ionizing front. 
Assuming that the main ionizing sources in the region 
are Objects la an d 2 (from Figure [6|), as sugges ted by 
Ortiz et al. (2007) and Roman-Lopes et al. ( 2009[ ), mag- 
netic field lines should be pointing toward these targets, 
which is indeed obtained from the polarization vectors 
mean orientation (which points mainly toward Objla). 

Such radial pattern is not noticed from the large-scale 
R-band map, although a large portion of this survey is 
composed by the ionized region, inside the Hll region. 
Furthermore, a perfect ring of polarimetric vectors is ob- 
viously not formed along the borders of the ionized area, 
as predicted by the simulations of the magnetic field lines. 
However, it is important to point out that the polarimet- 
ric survey is not a simple sky-projected slice of the Hii 
region. Rather, it represents the integrated effect of lin- 
ear polarization due to different interstellar structures 
superposed along the line-of-sight. For example, if we 
imagine the Hii region roughly as an expanding sphere, 
such expansion will be noticed spatially in the sky, but 
will also occur in the line-of-sight direction. The expan- 
sion of the outer borders will push and pull interstellar 
material in our direction, causing the observations of the 
ionized region to be composed by the sky-projected, in- 
tegrated polarizing effect due to the outer borders and 
the ionized region (if the observed star is located behind 
the Hii region). Therefore, the large-scale observations 
of the ionized region should account for this effect, and 
the radial pattern is not expected to be detected. 

However, even when considering such effects, an indi- 
cation of the dual behaviour of the polarimetric vectors' 
orientation inside and outside the ionized area is cer- 
tainly noticed when analysing the R-band and H-band 
polarization data. These observations give support to 
the results from the simulations, which provide similar 
qualitative predictions, at least when studying the spe- 
cific environment of the RCW41 region. 

5.2. The origin of the low polarization values within the 

Hll region 

In Section 14.31 we have seen that an overall lower R- 
band polarization degrees are detected inside the Hii re- 
gion, specially near the cluster. 
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Fig. 19. — {J ~ H) X {H — Ks) (color-color) diagrams using data 
from 2MASS, for two separate regions along RCW41 (as indicated 
by the red boxes): inside (a) and outside (&) the Hll region. This 
area is defined as being delimited by the radius of Rrise = 7.8it 1.0 
arcmin from the center, where a sharp rise in R- band polarization 
degree is detected, as previously shown in Figure [14] 



In order to distinguish between the two possible men- 
tioned causes responsible for such effect, two color-color 
diagrams have been built using 2MASS data, showing 
the distribution of stars inside (Figure 19 1) and outside 
(Figure 19 d) the limits defined by the sharp rise in polar- 
ization degree {Rrise = 7.8 ± 1.0 arcmin), which roughly 
coincides with the limits of the Hii region. A larger num- 
ber of stars is found outside the region, which is obviously 
due to the wider chosen area. However, when analyzing 
the distribution of stars along the reddening band (which 
gives an estimate of the interstellar extinction to each 
object), we notice that there are no major differences 
between both diagrams. In both areas, the majority of 
stars show a roughly uniform distribution of extinction 
values between Ay = and Ay ~ 12. Besides, we know 
that some concentrations mainly near the cluster (which 
comprises only a small area within the Hii region) may 
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pres ent larger mean local Ay values (as discussed in Sec- 
This analysis shows that no global decrease in 
liar extinction is detected inside the area. 

it seems that the observed decrease in po- 



a) R band Pol. Map; b= ( Ifi.l ± 0.4)": <B '>"7B = 0.202 ± 0.005 



3.3) 



tion 

interste . 
Therefore, 

larization degree inside the Hii region is probably caused 
by a lower polarization efhciency of the local interstellar 
medium. This may be explained by an intrinsic suppres- 
sion of the alignment of dust grains with respect to the 
interweaving magnetic field lines, which is expected given 
the dynamical and turbulent conditions of such interstel- 
lar environment. 

Furthermore, in the simulations performed by |Arthur| 
et al. (2011 ), it was shown that the intensity of magnetic 
heids in the photo- ionized area is lower {B = 0.1 — lO/iG), 
when compared with the neutral areas from the border 
of the Hii region {B = 10 — 30/iG). This is probably ex- 
plained due to the fact that the expansion of the ionized 
region causes the magnetic field lines to be pushed aside 
and therefore removed from inside the area, weakening 
the magnetic field in this location. If this is the scenario 
occurring in RCW41, a further depletion of polarization 
efficiency would be expected, contributing to the lower 
P values detected near the cluster. 

5.3. The turbulent component of the magnetic field 
toward RCW41 

The polarimetric dataset toward RCW41 is entirely 
suitable to apply a method which aims on estimating the 
turbulent contribution to the magnetic field li nes, rela- 
tive to its uniform component (Hildebrand et al.|2009p . It 
is based on the assumption that the interstellar magnetic 
field B(x) is composed of a large-scale structured field 
Bo(x), together with a turbulent component Bt(x). Ini- 
tially, the procedure consists basically in computing the 
second-order structure function of the polarization angles 
((A$^(Z)), hereafter SF), which is defined as the average 
value of the squared difference in polarization angles be- 
tween two p oints separated by a distance I (see Equation 
(5) given by [Falceta-Gongalves et al^|2008 ) . The square 
root of this quantity is denominated angular dispersion 
function (ADF), and provides information on the disper- 
sion of polarization angles as a function of the length 
s cale toward a specific inte rstellar area. 
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Hildebrand et al. (2009) pointed out that within the 
range d < I <^ d (where S and d are the correlation 
lengths respectively characterizing Bt(x) and Bo(x)), 
the SF is composed by: a constant turbulent contribu- 
tion (denoted by &); a smoothly increasing contribution 
with respect to the length I, as expected from the uni- 
form field (the slope of this linear behaviour is repre- 
sented by m); and an almost constant contribution due 
to measurement uncertainties (a{l)). Being statistically 
independent, these factors must be added quadratically, 
leading to 



(4) 



The ratio of the turbulent to large-scale magnetic field 
strength is related to the quantity b through the following 
expression: 
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Fig. 20. — Angular dispersion functions (ADF) constructed by 
using the polarization angles from the R band (a) and the H band 
(b) polarimetric datasets. Each bin (respectively of widths 2 and 
0.25 arcmin) have already been corrected due to the mean measure- 
ment uncertainties. The red dashed line denote the curve obtained 
through the linear fit of the SF, by using the first six points on each 
case (blue circles). The turbulent contribution b is determined by 
the interception of this curve at ( = 0, being subsequently used to 
compute the ratio of the turbulent to the large-scale magnetic field 
strength, which arc shown at the top of each diagram. 

Therefore, by constructing the SF (and computing the 
measurement uncertainties cr(^)) using the polarization 
angles, it is straightforward to apply a linear fit with 
Equation Hand determine the b factor (by the intercept 
at I = 0). 

Using the polarization angles from the R and H polari- 
metric datasets, we have computed the SF, as well as the 
mean measurement uncertainties cr(/). Thereafter, the 
ADF was built to each case, as shown in Figures [20^ and 
[20j 3, which are respectively related to the R and H band 
polarization data. Prior to calculating the ADF, the fac- 
tor cr'^{l) was subtracted within each bin of the SF, as 
suggested in Eq uation [4| (therefore , the value of each bin 
in Figure mis ^(A$2(7))_cr2(Z)). In dia grams 20a and 
[20j 3, bin widths of respectively 2 and 0.25 arcmin have 
been used. Only the first six points of the ADF (denoted 
by blue open circles) are used in the linear fit of Equa- 
tion |4] in order to guarantee the validity of the / <C d 
regime, and the fitted function is represented by the red 
dashed line. The ratio of the turbulent to large-scale 
magnetic field strength was computed through Equa- 
tion fs) resulting in {Bfy^^/Bo = (20.2 ± 0.5)% and 
{Bf)^^/Bo = (53 ± 1)% for the R and H band datasets, 
respectively. 

The calculation shows that the turbulent component of 
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the magnetic field is much stronger within the small scale 
area surrounding the cluster (i.e., where the H band data 
is distributed) , when compared with the global turbulent 
contribution considering the entire area surrounding the 
Hii region (the R band data distribution). In fact, a 
higher turbulence ambient is expected to exist nearer the 
cluster due to the larger interaction between the young 
stars and the interstellar medium within this area. We 
may also associate this evidence with our hypothesis that 
the low polarization values inside t he H ii region (mainly 
near the cluster, as seen in Figure 14 3), indeed occurs 
due to low polarization efhciency. The higher turbulence 
level probably represents the main factor inhibiting the 
alignment of dust grains with respect to the magnetic 
field lines, causing the lower polarization values. 

6. CONCLUSIONS 

Using two sets of NIR photometric data collected us- 
ing the SofI camera at the NTT 3.5 m telescope, we have 
characterized the stellar population of a deeply embed- 
ded stellar cluster related to the RCW41 Hii region. 
Furthermore, through optical (V, R and I bands) and 
NIR polarimetric observations (H band), obtained re- 
spectively with similar polarimetric modules at the 0.9 m 
telescope on CTIO/Chile and the 1.6 m telescope on 
OPD/Brazil, we have built linear polarization maps of 
the star-forming area. Our main results from both tech- 
niques are listed below: 

(i) Analysis of photometric colors of cluster's stars in 
comparison with field stars shows that a number of the 
cluster's objects present color excess, which is typically 
associated with circunstellar disks of young PMS stars. 
The cluster also presents several highly reddened objects 
(Ay > 20 mag); 

(ii) A mean visual extinction of Ay = 7.6 ± 2.0 is es- 
timated for the cluster's objects, based on the location 
of the majority of the stars on the color-color diagram, 
which roughly defines a band parallel to the reddened 
locus of CTT stars. This is corroborated by an indepen- 
dent Ay computation using MS stars from the cluster; 

(iii) The morphology of the cluster is divided between 
a main portion and a sub-cluster. There seems to exist a 
correlation between the location of the highly obscured 
objects with the HNCO contours, which is mostly con- 
centrated toward the north of the cluster's main part. 
Therefore, these are probably highly embedded objects, 
and a local rising extinction gradient in the south-north 
direction is probably present near the cluster. This ev- 
idence supports the idea that the interstellar material 
surrounding Obj2 (a 09v star at the sub-cluster, pre- 
sumably the most massive object in the field) have al- 
ready been swept out, while near the main portion, gas 
and dust are still beginning to be affected, according to 
the HNCO contours; 

(iv) The color-magnitude diagrams of the cluster's 
stars show that the majority of the objects are displaced 
toward higher color values, indicating that these stars are 
mainly young objects, also presenting color excess and 
non-uniform interstellar extinction. Comparison with 
PMS isochrones provide a mean cluster age between 2.5 
and 5.0 Myr. Stars from the turn-on region (more mas- 
sive) are mostly concentrated near the subcluster, sug- 
gesting that this area has been formed earlier, subse- 
quently triggering star formation at the main portion of 



the cluster. Results from item (iii) corroborate this hy- 
pothesis; 

(v) The nature of some objects previously observed 
through NIR spectroscopy was reviewed. Objlb, the 
close companion to the massive BOv central object, is 
highly reddened and presents spectroscopic features sim- 
ilar to T Tau Sb, being therefore a very young T Tauri 
star. Objects 3, 5, and 6 were previously identified as 
late-type field stars, although their position in the color- 
magnitude diagrams is consistent with the T Tauri clas- 
sification. Examples have been found of young objects 
with similar spectral features. Specifically, Obj3 is highly 
reddened {Ay > 20) and very luminous, with spectral 
features that may be consistent with a medium mass 
YSO. Its weak CO absorption lines may be explained 
due to veiling from the circunstellar dust emission; 

(vi) Using the density distribution of stars in the 
field, the cluster's center has been determined to 
be at (ac,^5c)20oo=(9''16'"43.5^-47°56'24.2"). A two- 
parameter King function has been adjusted to the radial 
density profile, considering both the entire area around 
the cluster and a second fit where the portion encom- 
passing the sub-cluster was removed, providing a much 
sharper profile. The general morphology shows that al- 
though some stellar scattering exist toward the south, the 
main cluster is probably of the centrally condensed-type. 
This shape is correlated to the spherical and similarly 
sized cloud core in which the cluster is immersed, as ob- 
served from the 1.2mm dust continuum emission toward 
RCW41. 

(vii) The R-band polarimetric survey of a large area 
encompassing the entire Hii region, was used to build im- 
ages of mean polarization degree and position angle. A 
radial analysis starting from the center of the ionized area 
revealed a sharp increase in polarization degree roughly 
coinciding with a conspicuous interstellar ring detected 
at the MSX 8.2^m emission, which delimits the ionized 
area. This evidence reflects a global decrease in polariza- 
tion degree inside the Hii region, mainly near the cluster, 
which is most probably due to low grain alignment effi- 
ciency caused by the turbulent environment and weak 
magnetic fields in this area; 

(viii) The image of mean polarization angle exhibit 
large areas, mainly toward the southwest of RCW41, 
where polarization angles are slanted outwards, follow- 
ing the borders of the ionized area. Such slanting coin- 
cides with the limits of the region indicated by the MSX 
8.2/im emission, and is supported by simulations of the 
expansion of Hii regions, which predicts that the mag- 
netic field lines at the neutral gas outside the limits of 
the region, tends to lie parallel to the ionization front; 

(ix) The H-band polarimetric survey, focused on a 
smaller area towards the cluster, shows that, although 
a large dispersion of polarization angles are present, 
a mean pattern may be noticed: polarization vectors 
roughly follow the direction of several interstellar fila- 
ments, emanating radially from the cluster, as revealed 
by the Spitzer 5.8/im image. Therefore, the mean orien- 
tation is roughly perpendicular to the ionization front, 
a result which is also in agreement with the predicted 
morphology of the magnetic field lines inside the ionized 
area, from simulations of expanding Hii regions; 

(x) By combining polarimetric optical data (V, R, and 
I) and the NIR H-band, the Serkowski function has been 
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fitted towards several stars mainly near the cluster. A 
large rotation of polarization angles is noticed toward 
some stars, revealing that the light beam has probably 
trespassed different interstellar layers along the line-of- 
sight. The total-to-selective extinction ratio has been 
computed toward each of the fitted objects, and a mean 
value of i? = 2.96 ±0.42 has been obtained, which agrees 
with the standard 3.09 value for the typical interstellar 
medium extinction law. 

(xi) Through the statistical correlation analysis (com- 
puting the ADFs using the R and H band polarimetric 
datasets), we have shown that the ratio of the turbulent 
to large-scale magnetic field strength is (20.2 ± 0.5)% 
when the entire area surrounding the Hll region is con- 
sidered, and (53 ± 1)% if only the area near the cluster is 
taken into account. The higher turbulence level within 
the cluster area reinforces the hypothesis that the low 
polarization values at this region occurs due to the lower 
polarization efficiency. 
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APPENDIX 

PHOTOMETRIC AND POLARIMETRIC DATA 

The lists related to the photometric data are available in Tables [I] (Dataset A) and |2j (Datsct B), and the complete 
multi-band polarimetric sample is provided in Table [3j Tables [l] [5] and [3] are publishedin its entirety in the electronic 
edition of the Astrophysical Journal. A portion of each table is shown here for guidance regarding its form and content. 



TABLE 1 

NIR PHTOMETRIC DATA OF RCW41 FROM DATASET A. 





'52ooo(° ) 


J 


AJ 


H 


AH 


Ks 


AKs 


9 16 47.94 


-47 57 17.8 


9.712 


0.022 


8.962 


0.015 


8.653 


0.026 


9 17 1.91 


-47 56 44.7 


10.999 


0.020 


9.456 


0.018 


8.786 


0.016 


9 16 43.49 


-47 56 22.9 


10.684 


0.021 


9.865 


0.015 


9.492 


0.026 


9 17 4.34 


-47 55 49.5 


12.801 


0.020 


10.852 


0.018 


9.959 


0.016 


9 16 47.67 


-47 57 21.2 


12.037 


0.022 


10.910 


0.015 


10.486 


0.026 


9 16 41.89 


-47 56 15.4 


14.555 


0.038 


12.041 


0.015 


10.627 


0.026 


9 16 47.05 


-47 56 44.0 


12.636 


0.021 


11.481 


0.015 


11.039 


0.026 


9 17 8.72 


-47 57 16.1 


12.982 


0.020 


11.764 


0.019 


11.183 


0.016 


9 16 45.57 


-47 56 56.0 


12.387 


0.021 


11.707 


0.015 


11.391 


0.026 


9 16 44.89 


-47 55 38.6 


13.740 


0.021 


12.407 


0.015 


11.483 


0.026 



Note. — The columns above respectively represent each of the objects' equatorial coordinates (a, (5), as well the .7, H and Ks magnitudes, 
together with each uncertainty (A J, AH and AKs). Values marked with * represent undetected stars (in each specific band), or saturated objects. 



TABLE 2 

NIR PHTOMETRIC DATA OF RCW41 FROM DATASET B. 



a20oo('""°) <520oo(°'") H AH Ks AKs 



9 16 47.95 


-47 57 17.8 






8.621 


0.036 


9 17 4.88 


-47 54 4.8 


10.411 


0.016 


9.534 


0.028 


9 16 43.49 


-47 56 22.8 


9.944 


0.024 


9.543 


0.028 


9 17 4.34 


-47 55 49.4 


10.870 


0.024 


9.995 


0.026 


9 16 47.68 


-47 57 21.2 


10.920 


0.016 


10.455 


0.021 


9 16 56.95 


-47 58 23.5 


11.194 


0.015 


10.500 


0.016 


9 16 41.88 


-47 56 15.4 


12.068 


0.022 


10.666 


0.024 


9 16 25.46 


-47 59 30.9 


11.337 


0.014 


10.765 


0.011 


9 16 34.57 


-47 56 58.3 


11.650 


0.015 


10.814 


0.015 


9 16 50.51 


-47 53 16.6 


11.093 


0.035 


10.924 


0.029 



Note. — The columns above respectively represent each of the objects' equatorial coordinates (a, (5), as well the H and Ks magnitudes, together 
with each uncertainty {AH and AKs). Values marked with * represent undetected stars (in each specific band), or saturated objects. 
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TABLE 3 

The V, R, I and H polarimetric data set. 





<^200o(° ) 


Band 


P{%) 






a0{°) 


9 15 37.95 


-48 1 24.9 


R 


2.70 


0.44 


144.6 


4.7 


9 15 41.72 


-47 53 10.7 


R 


5.33 


1.26 


3.3 


6.8 


9 15 43.66 


-47 47 39.1 


R 


10.62 


1.86 


1.8 


5.0 


9 16 22.63 


-47 57 57.1 


V 


10.47 


3.34 


171.1 


9.1 






R 


7.97 


1.44 


163.9 


5.2 






I 


6.12 


0.99 


162.1 


4.6 


9 16 25.00 


-47 58 1.0 


R 


10.72 


1.76 


177.8 


4.7 






I 


7.32 


0.75 


0.1 


3.0 


9 16 36.83 


-47 56 27.8 


H 


12.06 


1.59 


111.5 


6.3 


9 16 38.41 


-47 57 14.9 


V 


3.16 


0.31 


167.7 


2.9 






R 


2.51 


0.28 


168.5 


3.2 






I 


2.86 


0.35 


174.1 


3.5 






H 


0.87 


0.07 


177.8 


5.5 


9 16 40.99 


-47 55 9.9 


V 


3.60 


0.61 


138.7 


4.9 






R 


2.76 


0.16 


148.0 


1.7 






I 


2.35 


0.09 


150.8 


1.1 






H 


1.03 


0.07 


155.4 


5.3 


9 16 42.57 


-47 57 24.4 


V 


4.67 


5.42 


159.6 


33.2 






R 














I 


2.66 


2.53 


139.4 


27.2 






H 


3.99 


1.09 


165.4 


9.3 


9 16 44.28 


-47 58 17.1 


V 


0.85 


0.38 


178.0 


12.7 






R 


0.78 


0.37 


173.6 


13.7 






I 


0.44 


0.34 


170.7 


22.4 






H 






** 




9 16 45.47 


-48 3.4 


H 


11.43 


1.73 


147.2 


6.6 



Note. — The columns above respectively represent each of the objects' equatorial coordinates (a, 5), the spectral band of the polarization 
measurement (V, R, I, or H), the polarization degree and its uncertainty (P and crp), and the polarization angle with its uncertainty {0 and aO). 
Polarization and angle values marked with ** represent objects that were detected in the indicated band, but presents bad polarization values, as 
defined in Scction |2.2| 



